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Angular  momentum  is  prima  facie   a  formidable  obstacle  in  the  theory 
of  star  formation:  without  rotational  braking  during  star  formation, 
stars  would  rotate  with  speeds  close  to  that  of  light.  The  present 
investigation  suggests  that  magnetic  torques  acting  on  a  rotating, 
contracting,  cool  interstellar  cloud  which  is  permeated  by  a  frozen- in 
magnetic  field  coupling  the  cloud  to  its  surroundings,  rotationally 
decelerate  a  cloud,  constraining  it  to  co-rotate  with  the  background 
medium.  Angular-  momentum  is  thus  efficiently  transferred  from  a  collapsing 
cloud  to  its  surroundings. 

We  examine  angular  momentum  transfer  from .cool,  rotating,  stellar-mass 
condensations,  collapsing  isothermally  or-  a  magnetically-diluted  dynamic 
time  scale.  Some  mechanisms  are  discussed  for  forming  gravitationally- 
bound  protostellar  condensations  within  cool,  dense,  molecular  clouds. 
Rotation  induces  a  toroidal  magnetic  field  and  the  accompanying  magnetic 
stresses  generate  a  set  of  Alfven  waves  which  propagate  into  the  background 
medium,  thereby  transporting  angular  momentum  from  a  cloud  to  its  sur- 
roundings. A  modified  virial  approach  is  employed  to  calculate  time- 
dependent  quantities  of  interest  at  the  cloud's  surface  in  order  to 
estimate  the  braking  efficiency  of  the  magnetic  torques. 


It  is  found  that  so  long  as  a  cloud  remains  magnetically  coupled 
to  its  surroundings,  the  magnetic  torques  constrain  a  cloud  to  co-rotate 
with  the  background  medium.  Centrifugal  forces  are  always  kept  well 
below  gravity.  The  one  single  factor  most  important  in  determining 
the  angular  momentum  of  a  protostar  is  the  ionization  rate  in  dense 
magnetic  clouds:  the  degree  of  ionization  controls  the  coupling  of  a 
cloud  to  the  galactic  magnetic  field.  The  fractional  ionization  in 
dense  magnetic  clouds  is  therefore  discussed  in  some  detail. 

The  angular  momentum  of  magnetically-braked  protostars  is  shown  to 
be  consistent  with  the  observed  angular  momenta  of  close  binary  systems 
and  single  early-type  main-sequence  stars.  The  hypothesis  of  magnetic 
braking  offers  support  to  the  fission  theory  for  the  formation  of  close 
binary  systems,  and  is  able  to  account  for  the  relative  paucity  of  single 
stars.  The  calculations  also  suggest  a  common  mode  of  formation  for 
(close)  binary  and  planetary  systems. 

This  investigation  shows  that  magnetic  fields  do  indeed  play  an 
important,  if  not  dominant,  role  during  the  early  stages  of  star 
formation.  Detailed  numerical  hydrodynamic  collapse  models  have,  as  yet, 
ignored  the  possible  effects  that  magnetic  fields  may  have  on  the 
structure  and  evolution  of  a  protostar.  Such  models  are  therefore  highly 
suspect  and  probably  not  physically  realistic. 


SECTION  I 
INTRODUCTION 

Angular  Momentum  Problem 

Traditionally,  magnetic  fields  and  angular  momentum  have  presented 
formidable  problems  to  the  theory  of  star  formation  (cf.  Mestel  1965). 
Due  to  the  high  conductivity  of  the  interstellar  medium,  the  frictional 
coupling  between  plasma  and  neutral  gas  is  sufficient  to  cause  the  large- 
scale  galactic  magnetic  field  to  become  'frozen'  into  the  fluid  and 
dragged  along  with  it.  Accordingly,  the  magnetic  energy  density  of  a 
collapsing  interstellar  cloud  (or  fragment)  increases  as  the  cloud  contracts, 
arid  the  collapse  is  retarded  and  subsequent  fragmentation  may  be  prevented. 
Condensations  in  the  interstellar  medium  will  also  possess  angular 
momentum  by  virtue  of  local  turbulence  or  galactic  rotation.  A  simple 
calculation  shows  that  a  main- sequence  star  would  rotate  with  an  equa- 
torial speed  close  to  that  of  light  if  it  were  formed  by  isotropic 
compression  from  the  interstellar  gas,  conserving  angular  momentum  during 
contraction.  Of  course  it  is  doubtful  that  stars  could  ever  form  under 
such  conditions  since  centrifugal  forces  at  the  equator  will  increase 
faster  than  the  gravitational  forces,  ultimately  resulting  in  a 
rotational  instability. 

Mestel  and  Spitzer  (1956)  and  Nakano  and  Tademaru  (1972)  have 
shown  that  the  'magnetic  field  problem'  is  only  temporary.  Affibipolar 


■ 


diffusion  allows  the  field  to  uncouple  from  the  gas  when  the  fractional 
ionization  is  reduced.  Furthermore,- Mouschovias  (1976a, 1976b)  has  shown 
that,  at  least  for  relatively  low  gas  densities,  some  material  may  stream 
preferentially  down  the  magnetic  field  lines,  thereby  increasing  the 
ratio  of  gravitational  to  magnetic  energy  within  a  condensation. 

Radio- frequency  observations  of  molecular  clouds  do  not  show  any 
clouds  rotating  much  faster  than  the  Galaxy  (e.g.  Heiles  1970;  Heiles  and 
Katz  1976;  Bridle  and  Kesteven  1976;  Kutnex,  et  al .   1976;  Loren  1977,  and 
private  communication;  Lada,  et  al.   1974).  Main-sequence  stars  are 
observed  to  rotate  with  equatorial  velocities  ranging  from  a  few  hundred 
kilometers  per  second  for  the  early-type  stars  to  just  a  few  kilometers 
per  second  for  stars  later  than  spectral  type  F5  (Struve  1930;  Abt  and 
Hunter  1962).  Evidently,  nature  has  found  a  solution  to  the  angular 
momentum  problem. 

A  variety  of  mechanisms  have  been  proposed  to  reduce  the  angular 
momentum  of  collapsing  clouds  and  protostars.  Hoyle  (1945)  and  McCrea 
(1960,  1961)  have  suggested  that  condensation  may  take  place  in  regions 
where  the  local  turbulence  is  abnormally  small.  However,  each  object  is 
still  likely  to  have  somewhat  more  angular  momentum  than  is  found  in 
single  main-sequence  stars.  Preferential  mass  flow  along  the  rotational 
axis  would  increase  the  gas  density  at  constant  angular  velocity.  However, 
this  process  is  not  without  its  own  difficulties  (Mestel  1965;  Spitzer  1968a). 
More 'attention  has  been,  given  to  the  possibility  of  transforming  the 


Ambipolar  diffusion  ordinarily ' refers  to  the  process  of  charged- 
particle  diffusion  due  to  a  balance  between  a  space-charge  electric  field 
and  density  gradients  (of.  Krall  and  Trivelpiece  1973).  In  the  astro- 
physical  literature,  ambipolar  diffusion  refers  to  the  drift  of  a  weakly 
ionized  plasma  across  a  magnetic  field. 


intrinsic  (spin)  angular  momentum  of  a  single  massive  protostar  into  the 
orbital  angular  momentum  of  a  multiple  star  system  (Larson  1972a;  Black 
and  Bodenheimer  1976).  However,  as  Mouschovias  (1977)  points  out,  the 
angular  momentum  of  such  a  hypothetical  system  is  still  some  two  orders 
of  magnitude  greater  than  that  observed  for  the  long-period  (visual) 
binaries.  Dicke's  (1964)  claim  that  the  interior  of  the  Sun  is  in 
rapid  (differential)  rotation  suggests  that  single  stars  may  store  a 
large  amount  of  angular  momentum  beneath  their  surface.  Although  not 
accounting  for  the  possible  stabilizing  effect  of  toroidal  magnetic  fields, 
Goldreich  and  Schubert  (1967)  have  shown  that  a  necessary  condition  for 
stability  in  differentially  rotating  stars  of  homogeneous  chemical 
composition  is  that  the  specific  angular  momentum  (i.e.  angular  momentum 
per  unit  mass)  should  increase  with  increasing  distance  from  the  rotational 
axis.  Thus  it  appears  unlikely  that  a  differentially  rotating  main- 
sequence  star  can  have  an  angular  momentum  much  in  excess  of  a  uniformly 
rotating  star.  Furthermore,  convective  mixing  and  poloidal  magnetic 
fields  redistribute  angular  momentum  in  the  direction  of  rigid-body 
rotation. 

It  has  often  been  suggested  that  the  angular  momentum  of  a  con- 
tracting cloud  or  protostar  may  not  be  conserved.  That  is,  angular 
momentum  may  be  transferred  in  some  manner  to  the  surrounding  interstellar- 
material.  Weizsacker  (1947)  has  argued  that  a  rapidly  rotating  star  ;^ill 
be  rotationally  decelerated  as  angular  momentum  is  transferred  from  the 
star  to  its  surroundings  by  turbulent  viscosity.  Ter  Haar  (1949) 
subsequently  showed  that  Weizsacker 's  purely  hydrodynamic  mechanism  for 
angular  momentum  transport  is  probably  not  very  efficient.  Recently, 
Sakurai  (1976)  has  calculated  the  braking  torque  on  a  Jacobian  ellipsoid 


by  a  tidal  acoustic  wave  which  is  generated  in  the  surrounding  medium 
by  the  rotating  configuration.  However,  as  Sakurai  points  out,  the 
effectiveness  of  the  braking  for  pre-main-sequence  stars  is  uncertain 
because  the  braking  time  is  of  the  same  order  of  magnitude  as  the  time 
scale  of  evolution. 

Hydromagnetic  braking  appears  to  be  more  efficient  in  disposing  of 
angular  momentum.  In  an  attempt  to  account  for  the  sun's  observed  slow 
rotation,  Alfven  (1942)  suggested  that  the  interaction  of  the  sun's  dipole 
magnetic  field  with  the  surrounding  'ion  cloud'  would  produce  a  torque 
on  the  sun  tending  to  brake  its  rotation.  Ter  Haar  (1949)  generalized 
this  concept  to  include  all  stars  magnetically  coupled  to  HII  regions. 
Lust  and  Schliiter  (195S)  examined  in  some  detail,  particularly  for  the 
special  case  of  torque- free  magnetic  fields,  the  transport  of  angular 
momentum  by  magnetic  stresses  acting  on  a  rotating  star. 

Hoyle  (1960)  has  proposed  three  stages  of  development  for  star 
formation:  (1)  the  initial  stage  when  a  condensation  is  magnetically 
coupled  to  its  surroundings  by  the  frozen- in  galactic  magnetic  field. 
Angular  momentum  is  efficiently  transferred  from  the  contracting  conden- 
sation to  the  surrounding  medium  with  the  condensation  being  constrained  to 
co-rotate  with  the  surroundings;  (2)  a  subsequent  phase  when  the  fractional 
ionization  becomes  low  enough  so  that  the  condensation  uncouples  from  the 
galactic  magnetic  field  via  ambipolar  diffusion,  afterwhich  angular  momentum 
is  effectively  conserved;  and  (3)  a  recoupling  with  the  galactic  field 
during  the  final  stage  of  slow  contraction  to  the  main  sequence.  Hoyle 
was  able  to  explain  the  anomolous  distribution  of  angular  momentum  within 
the  solar  system  (98%  of  the  total  angular  momentum  of  the  solar  system 
is  concentrated  in  the  planets  which  comprise  less  than  It  of  the  total 


mass  of  the  system)  as  being  the  result  of  a  hydromagnetic  transfer  of 
angular  momentum  from  the  primitive  solar  nebula  to  the  planetary 
material.  Hoyle's  calculations  were  confirmed  in  a  more  quantitative 
fashion  by  Dallaporta  and  Secco  (1975) . 

Following  Hoyle's  (1960)  paper  on  the  origin  of  the  solar  system 
(for  a  review  of  this  and  other  theories  of  solar  system  formation,  see 
Williams  and  Cremin  1968),  it  was  generally  believed  (McNally  1965; 
Huang  1973,  and  references  cited  therein)  that  single  main-sequence  stars 
of  spectral  type  F5  and  later  were  likely  to  have  planetary  systems,  and 
that  their  observed  slow  rotation  was  thus  explained  ipso  facto. 
Schatzman  (1962)  pointed  out  that  the  transition  between  stars  with  deep 
envelopes  in  radiative  equilibrium  and  those  with  well-developed  sub- 
surface hydrogen  convection  zones  occurred  among  the  F  types.  He 
introduced  a  theory  in  which  the  gas  emitted  by  jets  and  flares  associated 
with  the  active  chromospheres  of  the  later-type  stars  (those  stars  having 
subphotospheric  convective  zones)  is  magnetically  constrained  to  co-rotate 
with  the  star  out  to  very  large  distances  where  it  carries  away  a  large 
amount  of  angular  momentum  per  unit  mass.  This  theory  is  consistent  with 
observational  evidence.  T  Tauri  stars  undergoing  pre-main-sequence 
contraction  are  ejecting  matter  (Herbig  1962;  Kuhi  1964,  1966;  see, 
however,  Ulrich  1976).  The  observations  of  Wilson  (1966)  and  Kraft  (1967) 
show  a  connection  between  the  rotation  of  stars  and  their  age  as 
determined  by  chromospheric  activity  (measured  by  the  presence  of  H  and 
K  emission  lines  of  Call)  which  is  usually  associated  with  the  hydrogen 
convection  zone. 

Dicke  (1964) ,  Brandt  (1966) ,  Modisette  (1967) ,  and  Weber  and  Davis 
(1967)  have  calculated  the  solar-wind  induced  torque  on  the  Sun.  They 


conclude  that  the  torque  is  sufficient  to  halve  the  sun's  rotation  on 
a  cosmological  time  scale.  Elaborating  on  the  ideas  of  Schatzman, 
Mestel  (1968)  lias  formulated  a  theory  of  magnetic  braking  by  a  stellar 
wind.  Using  Mestel 's  results,  Schwartz  and  Schubert  (1969)  have  shown 
that  the  Sun  may  liave  lost  a  considerable  amount  of  angular  momentum 
if  it  passed  through  an  active  T  Tauri  stage.  Assuming  that  stars  in 
the  pre-main- sequence  stage  are  wholly  convective  (Hayashi  1961) , 
Okamoto  (1969,  1970)  has  shown  that  solar-type  stars  may  lose  almost  all 
of  fheir  angular  momentum  via  a  Schatzman-type  braking  mechanism  during 
pre-main-sequence  contraction. 

The  Schatzman-type  magnetic  braking  mechanism  would  apply  only  to 
those  stars  having  appreciable  subsurface  convection  zones  and  therefore 
enhanced  surface  activity  (e.g.  mass  loss).  This  may  account  for  the 
break  in  stellar  rotation  on  the  main  sequence  at  spectral  type  FS, 
although  as  mentioned  earlier,  it  may  be  that  in  some  cases,  angular 
momentum  has  been  transferred  to  a  surrounding  planetary  system.  It  is 
not  clear  that  early-type  stars  ever  develop  a  fully  convective  structure 
during  pre-main-sequence  contraction  (Larson  1969,  1972b).  Accordingly, 
for  these  stars  in  particular,  we  must  examine  the  possibility  of 
rotational  braking  during  the  early  pre-opaque  stages  of  star  formation. 

Ebert,  et  al .  (1960;  see  also  Spitzer  1968b  and  Rose  1973),  in  a 
pioneering  study  have  investigated  the  transfer  of  angular  momentum 
from  a  contracting  interstellar  cloud  which  is  magnetically  linked  to  the 
surrounding  interstellar  medium  by  the  frozen- in  galactic  magnetic  field. 
Kinks  in  the  field  lines  introduced  by  the  rotation  of  the  cloud  propagate 
into  the  surrounding  medium  in  the  form  of  magnetohydrodynamic  (MHD)  waves 
(in  this  case,  the  transverse  Alfven  mode  is  excited),  thereby  rotationally 
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decelerating  the  cloud.  This  mechanism  is  expected  to  be  operative  only 
so  long  as  the  cloud  remains  magnetically  coupled  to  the  background. 
As  the  collapse  proceeds  to  higher  densities,  ambipolar  diffusion 
(Mestel  and  Spitzer  1956;  Nakano  and  Tademaru  1972),  MHD  instabilities 
(Mestal  1965) ,  or  perhaps  intense  Ohmic  dissipation  (Mestel  and 
Strittmatter  1967)  may  act  to  uncouple  the  cloud's  field  from  the 
surrounding  medium.  Although  their  results  remain  somewhat  tentative  due 
to  the  uncertainties  in  the  formulation  of  the  problem  (e.g.  assumed 
cylindrical  symmetry) ,  it  appears  that  the  magnetic  torques  may  be 
sufficient  to  brake  the  cloud's  rotation  so  that  Hoyle's  (1960)  argument 
for  efficient  angular  momentum  transfer  during  the  initial  stage  (Hoyle's 
stage  1)  of  star  formation  is  supported.  In  a  more  detailed  general 
analysis,  Gillis  et  al.   (1974)  find,  in  one  particular  application  of 
their  somewhat  artificial  time -independent  pseudo-problem,  that  the 
magnetic  braking  is  "embarrassingly  efficient"  although  they  admit  that 
their  mathematical  approximations  introduce  some  degree  of  uncertainty. 

Kulsrud  (1971)  has  calculated  the  rate  of  emission  of  energy  in  the 
form  of  MHD  waves  (specifically,  the  fast  magnetosonic  mode)  for  a 
rotating,  time-dependent,  point  magnetic  dipole.  Kulsrud  has  shown  that 
stars  with  very  large  magnetic  fields  (e.g.  magnetic  A  stars)  and 
initially  small  rotation  may  be  decelerated  to  very  long  periods.  Indeed, 
this  mechanism  may  explain  the  anticorrelation  of  rotational  velocity  and 
surface  magnetic  field  strength  observed  for  the  magnetic  stars 
(Landstreet  et  al.   1975;  Hartoog  1977).  The  magnetic  accretion  theory 
of  Havnes  and  Conti  (1971)  and  the  centrifugal  wind  theory  of  Strittmatter 
and  Norris  (1971)  have  also  been  proposed  to  account  for  the  long-period 
Ap  stars.  Nakano  and  Tademaru  (1972),  Fleck  (1974),  and  Fleck  and 


Hunter  (1976)  have  adapted  Kulsrud's  result  (even  though  Kulsrud's 
formulae  are  strictly  applicable  only  to  a  periodically  time-varying 
dipolar  field)  in  order  to  estimate  the  efficiency  of  braking  for 
collapsing  interstellar  clouds.  The  results  of  Fleck  and  Hunter  are  in 
good  agreement  with  observations  of  molecular  clouds  and  stellar  rotation 
on  the  main  sequence. 

Prentice  and  ter  Haar  (1971;  see  also  Krautschneider  1977)  have 
suggested  that  a  collapsing  grain-cloud  may  lose  angular  momentum  to 
the  neutral  gas  component.  However,  this  mechanism  assumes  that  the 
grains  are  electrostatically  neutral,  and  it  ultimately  relies  on  a  hydro- 
magnetic  transfer  of  angular  momentum  to  the  outside. 

Present  Work 
The  purpose  of  the  present  investigation  is  to  show  that  magnetic 
fields  do  indeed  play  an  important ,  if  not  dominant ,  role  during  the 
early  stages  of  star  formation.  We  examine  angular  momentum  transfer 
from  a  cool,  rotating,  magnetic  cloud,  magnetically  coupled  to  its 
surroundings  prior  to  the  epoch  of  ambipolar  diffusion,  and  undergoing 
essentially  pressure-free  collapse  on  a  magnetically-diluted  dynamic 
time  scale.  Rotation  induces  a  toroidal  magnetic  field  in  the  neighborhood 
of  the  cloud  and  the  accompanying  magnetic  stresses  produce  a  net  torque 
acting  on  the  cloud  tending  to  keep  the  cloud  in  a  state  of  co-rotation 
with  its  surroundings.  We  do  not  attempt  a  detailed  solution  of  the 
coupled  hydrodynamic  and  electrodynamic  equations  as  to  do  so  would 
require  a  sophisticated  computer  code  to  handle  the  problem  numerically. 
Such  a  formidable  (if  not  impossible)  task  is  hardly  justifiable  in 
view  of  our  lack  of  understanding  of  many  of  the  details  of  the  star 
formation  process.  Instead,  we  employ  a  modified  virial  approach  to 
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calculate  time-dependent  quantities  of  interest  at  the  surface  of  a 
cloud  in  order  to  estimate  the  efficiency  of  the  magnetic  torques  in 
de-spinning  the  cloud.  Ke  compare  our  results  with  obsei-ved  properties 
of  molecular  clouds,  the  specific  angular  momenta  of  (close)  binary 
systems,  the  angular  momentum  of  the  protosun,  and  with  stellar  rotation 
on  the  main  sequence. 

Uncertainties  in  some  of  the  physical  processes  of  star  formation  and 
complexities  in  the  mathematical  formulation  of  the  problem  do,  of  course, 
necessitate  some  degree  of  approximation  and  simplification  in  order  that 
the  problem  remain  tractable.  We  cannot  hope  to  improve  on  the  approximate 
nature  of  any  theoretical  study  of  star  formation  until  we  better 
understand  the  observations  that  are  just  now  becoming  available. 


SECTION  II 
MAGNETIC  BRAKING 

Magnetic  and  Velocity  Fields 
It  has  been  established  (Heiles  1976,  and  references  cited  therein) 
that  a  lai-ge- scale  magnetic  field  pervades  the  T3alaxy.  Due  to  the  high 
conductivity  of  the  interstellar  medium,  this  field  is  'frozen'  into 
the  fluid  (Mestel  and  Spitzer  1956).  Consider  a  uniform,  spherical, 
interstellar  cloud  with  radius  R  which  is  contracting  isotropically.* 
Strict  flux  conservation  implies  that  the  magnetic  field  strength  B 
within  a  radially  contracting  cloud  increases  according  to 

B  =  Bo(Ro/R)2   ,  (1) 

where  the  subscripts  denote  initial  values.  As  a  consequence  of  flux- 
freezing  during  an  isotropic  collapse,  the  initially  uniform  (galactic) 
field  lines  are  drawn  out  from  the  cloud  into  a  nearly  radial  structure 
(Mestel  1966).  Accordingly,  we  approximate  the  magnetic  field  outside 
the  cloud  by  the  spherical  polar  coordinates  B  =  (B.,B„,B,)  where 

Rn 
Br  =  V~T  +1)cos9  (Z) 


Observations  of  condensations  in  the  interstellar  medium  spanning  a 
range  in  mass  from  the  massive  molecular  cloud  complexes  down  to  the  stellar- 
mass  Bok  globules  (Zuckerman  and  Palmer  1974,  and  references  cited  therein) 
indicate  an  approximate  spherical  geometry.  Isotropic  contraction  will 
be  partly  justified  and  partly  relaxed  in  a  later  section  of  this  paper. 
Of  course,  the  simplifying  assumption  that  a  cloud  is  uniform  and  has  a 
well-defined  boundary  at  R  is  somewhat  artificial  although  it  does  simplify 
the  calculations  and  it  is  not  expected  to  affect  the  validity  of  the  results. 
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B„  =  -B  sin9  (3) 

0    o 

B -  0   .  (4) 

For  r»R  ,  the  field  becomes  uniform  and  is  described  by  the  equation 
o 

t    =  B  (cos9,-sin6,0)  .  (5) 

0     0^     ' 

The  velocity  fields  outside  a  radially  contracting,  rotating  cloud 
are  given  by 

v  «  (vr,0,v^  ,  (6) 

where 

vr  -  R(r/R)n  (7) 

and 

v.  *  uirsinfl  .  W 

9 

In  Eq.  (?),!  =  $,  the  collapse  velocity  at  the  cloud  surface  (r=R) , 
and  we  set  the  exponent  n=l  in  accordance  with  the  findings  of  Gerola 
and  Sofia  (1975)  and  Fallon  et  al.  (1977,  and  private  communication) 
for  the  Orion  A  molecular  cloud.  However,  the  exact  value  of  the 
exponent  is  somewhat  uncertain  (of.  Loren  et  al.   1973;  Loren  1975,  1977; 
Snell  and  Loren  1977).  In  Eq.  (8),  u=u(r,t)£  is  the  angular  velocity  of 
the  material,  and  we  have  taken  the  axis  of  rotation  to  be  parallel  to 
I  ,  i.e.,  io=B  =£,  the  unit  vector  along  the  positive  z-axis.  The  z-axis 
thus  becomes  the  axis  of  symmetry  and  derivatives  with  respect  to  the 

azimuthal  coordinate  vanish,  i.e.,  gr  =  °  ■  Mouschovias  (private 
communication)  is  investigating  magnetic  braking  for  the  case  of  wi3Q, 
and  believes  that  the  braking  may  be  more  efficient  in  this  case. 
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Paris  (1971)  has  shown  a  tendency  for  the  torques  exerted  by  an  undetached 
magnetic  field  to  rotate  the  angular  momentum  vector  into  parallelism 
with  the  overall  direction  of  the  field.  Thus,  our  assumption  that 
u-E  =1  probably  more  closely  approximates  reality.  Velocities  at  the 
cloud's  surface  can  be  found  by  setting  r=R.  The  cloud  is  assumed  to 
rotate  rigidly  at  a  uniform  rate  u(riR)=u(R) .  The  magnetic  viscosity 
due  to  the  cloud's  frozen- in  magnetic  field  constrains  the  cloud  to  rotate 
uniformly  as  long  as  the  travel  time  of  an  Alfven  wave  through  the  cloud 
is  less  than  the  collapse  time. 

Torque  Equation 
The  shear  at  R  due  to  the  cloud's  rotation  generates  a  toroidal 
field  B  ,  and  the  resulting  magnetic  torques  react  on  the  rotation  field. 
The  magnetic  stresses  acting  to  minimize  B  generate  a  set  of  Alfven  waves 
which  propagate  into  the  surrounding  medium,  thereby  transporting  angular 
momentum  from  the  cloud  to  its  surroundings.  As  pointed  out  by  Liist 
and  Schliiter  (1955;  see  also  Mestel  1959),  the  magnetic  torque  exerted 
on  currents  within  a  given  volume  can  be  described  by  means  of  a  tensor 


\i =  EkifiV 


analogous  to  the  Maxwell  stress  tensor* 


(9) 


M 


'*■•> 


\» -%*»--*&  ■  (10) 

where  e,-.  and  6,  are,  respectively,  the  Levi-Civita  tensor  and 
Kronecker  delta,  and  x  ,  a=i,j,k  ,  denotes  the  Cartesian  coordinates. 


We  employ  the  Gaussian  system  of  electromagnetic. units. 
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If  (k,i,j)  is  a  cyclic  permutation  of  Eqs.  (9)  and  (10),  then 

\^*iVVi*  =  (xiVVi^-(xiyyi^  ■     cn) 

The  k-component  of  the  magnetic  torque  density  about  the  origin  is 
given  by 

so  that  the  magnetic  torque  acting  on  a  volume  V  may  be  transformed  into 
a  surface  integral: 

K*  ■  -  fDkAds  (13) 

where  n.  is  the  unit  normal  outward  from  the  surface  element  dS.  If  the 
surface  S  is  a  sphere  centered  on  the  origin,  then  the  total  outflow  of 
angular  momentum  is 

JV¥iSnM(-BrH*iVJVdS     '  (14} 

2 
where  B  is  the  radial  component:  although  the  magnetic  pressure  (B  /8tt) 

can  interchange  angular  momentum  between  field  streamlines,  only  the 

tension  along  the  field. lines  B.B./4tt  contributes  to  the  flux  across 

the  sphere  S  because  the  pressure  acting  normally  to  each  surface  element 

has  zero  moment  about  the  center.  The  total  torque  can  have  only  a 

z-component  since  our  system  is  symmetric  about  the  z-axis.  Letting  B 

denote  the  poloidal  component  of  the  magnetic  field  (i.e.,  B  =B  +Bg), 

the  total  torque  x   is  (Liist  and  SchlUter  1955) 

T  =  h  JVprsin8ds  -    ■  (15) 
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where  the  surface  element  for  a  sphere  of  radius  r  is  just 

dS  =  r2sin6d@d(|>  .  (16) 

For  a  rotation  field  described  by  Eq.  (8) ,  it  is  intuitively  clear 
that  the  toroidal  magnetic  field  vanishes  along  the  z-axis  where  v  is 
zero,  and  in  the  xy-p]ane  where  B  changes  sign.  Thus  one  can  write 
the  toroidal  field  as 

B  (r,6)  =  B.(r)sin9cos0  .  (17) 

Since  the  poloidal  magnetic  field  outside  a  collapsing  cloud  has  an 
almost  purely  radial  structure,  we  set  B  =B  so  that  using  Eqs.  (16)  and 
(17)  in  Eq.  (15)  and  carrying  out  the  appropriate  integration,  the  torque 
at  the  surface  of  the  cloud  becomes 

t  =  fa  R3B+(R)BSCR)  (18) 

where  B  (R)  is  the  surface  poloidal  field  and  is  given  by  Eq.  (1). 

Toroidal  Magnetic  Field 
The  time  dependence  of  a  frozen- in  magnetic  field  is  given  by 
(cf.   Jackson  1975) 

||=Vx(vxS)  .  (19) 

ox. 

Provided' that  the  ratio  |B,/B  |  does  not  become  large,  the  temporal 
behavior  of  the  surface  poloidal  field  should  be  adequately  described 
by  Eq.  (1)  if  R=R(t)  is  known.  For  the  v  and  B  fields  described  by 
Eqs.  (2),  (3),  and  (6) -(8),  the  time  dependence  of  the  toroidal  field, 
given  by  Eq.  (19)  is 
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^^Wv-rV^W1    '  (20) 


which  becomes 


J-ivv^ff^  -  «» 


where  we  have  made  use  of  the  convective  derivative, 
3t   at       8t    r  3r 


(22) 


In  Eq.  (21) ,  the  first  term  represents  the  convection  of  B,  due  to  v 
while  the  second  term  shows  clearly  the  expected  dependence  of  dB /dt 
on  the  shear  in  the  azimuthal  velocity  field  3co/3r  .  For  lvrl<<vA  > 
where 

v  -  -J-  (23) 

A  ^p 

is  the  Alfven  speed  in  a  plasma  having  a  mass  density  p,  the  convection 
term  is  unimportant  and  the  rate  of  growth  of  B,  is  determined  soley  by 
the  rotational  shear  3i,y'3r  . 

We  now  derive  an  approximate  expression  for  u(r) ,  and  finally, 
3u>/3r  .  The  equation  of  motion  in  a  fixed  non-rotating  inertial  frame  is 

P  !]      -  -V(Pg+P*)  ♦  1M  .  (24) 

'fixed 

Here,  P  is  the  thermal  gas  pressure,  *  denotes  the  gravitational 
potential  and 

V!**  (25) 

is  the  magnetic  force  density,  j  being  the  electric  current  density, 
and  c  is  the  speed  of  light.  The  transformation  of  v  between  a  fixed 
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frame  and  a  frame  rotating  with  angular  Velocity  to  is  given  by 
(cf.  Marion  1970) 


dv 


dv 
dt 


♦   0)XV 


(26) 


fixed  rotating 

Sue,  in  the  reference  frame  of  our  rotating  cloud  (u  being  the  cloud's 
angular  velocity)  the  equation  of  motion  for  the  velocity  field  given 
by  Eqs.    (6) -(8)  reads 


dv 

■at 


rotating 


-v'CP  +p*)+fM-  ^u-r-puxCuxr) 


(27) 


The  6 -component  of  this  equation  is 
P1IF=-VYP$)+V£ce    • 


(28) 


where  Vfi=S/38  ,  fMfl  is  the  0 -component  of  the  magnetic  force  density  and 


.-*•     „->   -K 


fc8  =  P[oo*(u*r)]e 

2 
=  pu  rsin6cos9 


(29) 


is  the  0-component  of  the  centrifugal  force  density.  The  first  term  on 
the  right-hand-side  of  Eq.  (28)  vanishes  for  a  spherically  symmetric  cloud. 
The  left-hand-side  of  the  equation  vanishes  as  well  since  the  velocity 
field  given  by  Eq.  (6)  assumes  v  =0  .  Thus,  Eq.  (28)  reduces  to 


^ffi   £c0  " 
Using  Ampere's  law 


(30) 


(31) 
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in  Eq.  (27) ,  the  magnetic  force  density  becomes 

£--&*«**■  (32) 

whence  the  6. -component 

B   ,      3B     B,    , 


Canbining  Eqs.  (29)  and  (33)  in  accordance  with  Eq.  (30),  using  Eqs.  (2) 
and  (3)  for  B  and  Bfl,  respectively,  yields  the  following  expression  for 


2 
a)  : 

R2 
w2  »  -^  [B2R2(-|  +l)+2B2(r)(sin2e-cos26)]  ,  (34) 

pr   °   r 

where  we  have  used  Eq.  (17)  to  write  out  the  explicit  r  and  6  dependence 
of  B,  .  The  effect  of  the  second  term  in  brackets  is  to  increase  oi  in 
the  equatorial  zones  (i.e.  the  xy-plane)  and  decrease  u  near  the  poles 
(i.e.  along  the  z-axis).  Averaged  over  a  sphere  of  radius  r  which  is 
concentric  with  the  cloud,  this  term  vanishes,  i.e., 


.•_2„  _  2,.  _  0 


/"(sin29-cos2e)d6 


<sin  6 -cos 


sV  =  J =  o  ,  (35) 


{« 


so  that  an  approximate  (average)  angular  velocity  for  the  material 
surrounding  the  cloud  is 

B  R   R2   , 

»  «  -§•£(-§  +D^  •  (36) 

r  p2  r 

This  procedure,  which  is  equivalent  to  neglecting  currents  in  the  radial 
direction,  is  similar  to  that  employed  by  Alfven  1967;  also  Alfveri  and 
Arrhenius  (1976)  in  deriving  his  'law  of  partial  co-rotation'  for  a 
magnetized  plasma. 
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The  gas  density  p  outside  the  cloud  will  be,  in  general,  some 
function  of  r.  The  theoretical  collapse  models  of  Hunter  (1969)  and 
the  observational  findings  of  Loren  (1977)  for  the  Mon  R2  molecular 
cloud  suggest 

p  =  Ps(r2  « 

where  p  =  3m/47tR3  is  the  density  at  the  surface  of  a  uniform  spherical 
cloud  having  a  mass  m.  Using  this  result  in  Eq.  (36)  and  differentiating 
with  respect  to  r  yields 


3m  _  ,4itR^  o  o 
3r  m{rSr       r2 


7  7        2 

B  R  r  R    i  R  K         v 
r      r   r 


(38) 


As  expected,  |^<0  •  Using  this  expression  for  3n>/3r  in  Eq.  (21), 
2 

K 

taking  -S»  1  (which  should  be  true  as  the  collapse  proceeds  and  has 

r 
the  virtue  of  somewhat  underestimating   the  rate  of  growth  of  B  initially) , 

and  evaluating  the  result  at  the  surface  of  the  cloud,  r=R,  yields 

^^=-fB^O)-(^^sinecose.  (39) 

This  is  a  linear  first-order  differential  equation  which  can  be  cast 
into  the  form 

An  integrating  factor  is  R  .  Making  use  of  the  free-fall  collapse 
velocity  (in  a  later  section,  we  modify  this  equation  to  take  into 
account  pressure  gradients  within  the  cloud) 

§'--[»<$-•&]*   ,  (41) 
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G  being  the  Newtonian  gravitational  constant,  and  defining 


n  =  R/Rn  , 


(42) 


the  solution  to  Eq.  (40)  assumes  the  form 


B.(R,6)  =  -n 


„  ,  B2R  sin8cos9  m 


111 


(1-1)- 


.   (43) 


Evaluating  the  integral  and  writing  B^(R,6)  -  B(jfR)sin9cose   in  accordance 
with  Eq.   (19)  gives 

1.4SX10"35B2R2       „     -A  7         35     2j  35     3, 


B.(R)  -  -n  B   (R  ) 

*    °        G  (m/mo) 

•l-(l-n)^" 


4n 


2    35    . 
""    T28£n 


.l+(l-n) 


(44) 


where  n    -  2.0*1033g  is  the  mass  of  the  Sun.     Asymptotically,  as   n-J), 
a 

BfR)+  r,"2   • 

Rotational  Deceleration 
The  net  torque  t  acting  on  a  rotating  cloud  is  related  to  the  time- 
rate-of-change  of  angular  momentum  by 


where 


dT 
dt  ' 


J  =  KmR  u 


(45) 


(46) 


is  the  cloud's  total  angular  momentum,  k  being  the  gyration  constant 
(k=0.4  for  a  homogeneous  uniformly  rotating  spherical  cloud).  Combining 
these  two  equations  yields  an  expression  for  the  rotational  deceleration 
of  the  cloud: 


l\) 


K    kMR 

The  collapse  velocity  at  the  cloud  surface  R  is  given  by  Eq.  (41),  and 
the  magnetic  torque  acting  on  the  cloud  is  determined  from  Eq.  (18) 
using  Eq.  (44)  to  evaluate  B (R) .  Notice  that  for  t=0,  the  above 
expression  reduces  to  angular  momentum  conservation.  Angular  momentum 
is  transferred  from  the  cloud  to  its  surrounding's  so  that  t  is  intrin- 
sically negative  and  the  cloud  is  rotationally  decelerated. 

A  discussion  of  some  of  the  relevant  time  scales  is  in  order.  From 
Eq.  (47)  it  is  apparent  that  braking  will  be  efficient  only  if  the  second 
term  on  the  right-hand-side  dominates  the  first.  Using  Eq.  (23)  for  the 
Alfven  speed,  one  can  easily  show  that  this  is  equivalent  to  the  following 
condition: 

\>   |vRv^|  ,  (48) 

where  vR=R  and  v,=u)R  is  the  cloud's  surface  rotational  velocity  in  the 

equatorial  zones.  The  combined  radial  and  azimuthal  mass  motion  must  not 

exceed  the  wave  speed  at  the  surface  if  the  magnetic  stresses  are  to 

transport  angular  momentum  to  the  surrounding  medium.  A  crude  estimate 

I 
of  the  power  radiated  away  via  MHD  waves  is  given  by 

E  ' 
P  ~     rot 

h 
=  |kMRvaoj2  ,  (49) 


where 


^ot^^2"2  CSO) 


is  the  cloud's  rotational  kinetic  energy  and 
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tA  -  R/vA  (51) 

is  a  measure  of  the  characteristic  hydromagnetic  time  scale,  i.e.,  the 
travel  time  for  an  Alfven  wave  traversing  the  cloud.  (Interestingly 
enough,  this  order- of -magnitude  estimate  for  the  power-loss  is,  excepting 
for  a  constant  of  order  one,  just  that  predicted  by  the  magnetic  braking 
model  of  Ebert,  et  al .    (I960)).  Since  P  =  -ra,  Eq.  (47)  becomes 

£-.-$.$♦«  .  (52) 

whence  the  condition 

vA  >  |vR|  (S3) 

in  order  that  braking  be  efficient.  A  measure  of  the  characteristic  time 
scale  for  free-fall  collapse  is 

t,  «  —  (54) 

f   VR 

so  that  the  condition  for  efficient  braking  becomes 

tA  <  tf  (55) 

For  a  marginally  unstable  cloud  collapsing  from  rest,  this  condition  is 
satisfied  during  the  initial  contraction  stage  since  v.  is  typically  a 
few  kilometers  per  second  in  the  interstellar  medium.  In  fact,  Houschovias 
(private  communication)  believes  that  the  magnetic  stresses  acting  on  the 
surface  of  a  contracting  cloud  will  prevent  vR  from  ever  exceeding  v. . 
It  is  sometimes  argued  that  because  the  magnetic  energy  of  a  gravitationally- 
bound  condensation  can  never  exceed  the  gravitational  energy,  the  travel 
time  of  Alfven  waves  througli  the  condensation  is  at  least  equal  to,  and 
may  well  be  considerably  longer  than,  the  free-fall  time  which  is  given  by 
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tf  =  ffy*    .  (56) 

However,  this  is  the  time  required  for  complete   collapse  to  a  zero- 
radius  singularity  (cf.   Hunter,  1962) .  It  is  more  appropriate  to  compare 
time  scales  of  interest  with  the  'instantaneous'  dynamic  time  scale  as 
given  by  Eq.  (54).  As  pointed  out  by  Mestel  (1965)  and  Mouschovias 
(1976a,  1976b,  1977) ,  the  free-fall  time  as  defined  by  Eq.  (56)  may  have 
but  an  academic  significance  for  clouds  with  frozen- in  magnetic  fields. 


SECTION  III 
STAR  FORMATION 

Shock- induced  Star  Formation 

The  fact  that  young  stars  are  frequently  found  in  clusters  suggests 

that  stars  are  formed  by  a  fragmentation  process  which  occurs  during 

the  gravitational  collapse  of  large  interstellar  clouds.  According  to 

the  Jeans  (1928)  instability  criterion,  the  minimum  tmstable  mass  is 

_3 
related  to  the  gas  temperature  T  (K)  and  particle  density  n  (cm  ) 

through  the  relation 

3J  >  10  (lV/2  •  <S7) 

in      vn 
e 

Due  to  the  isothermal  behavior  of  the  interstellar  medium  at  relatively 
low  gas  densities  (Gaustad  1963;  Gould  1964;  Hayashi  and  Nakano  1965; 
Hattori  et  al.   1969;  see  also  Appendix  A  of  this  paper),  the  minimum 
unstable  mass  decreases  as  the  collapse  proceeds  to  higher  gas  densities 
so  that  a  large  collapsing  cloud  is  expected  to  fragment  into  a  number  of 
smaller  stellar-mass  condensations. 

However,  rather  compelling  theoretical  arguments  and  observational 
evidence  have  been  presented  suggesting  that  gravitationally-bound 
stellar-mass  condensations  (i.e.  protostars)  may  form  directly  out  of 
the  interstellar  medium  without  recourse  to  fragmentation.  Ebert  (1955) 
and  McCrea  (1957;  see  also  the  discussion  in  Mestel  1965)  have  shown  that 
external  pressures  of  the  order  10  to  10  cm  "  K  can  reduce  the  minimum 
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unstable  mass  to  stellar  order.  Such  extreme  pressure  variations  are 
known  to  exist  in  the  interstellar  medium  (Jura  1975) . 

Shock  waves  propagating  in  the  interstellar  medium  can  increase 
the  gas  density  up  to  two  orders  of  magnitude,  and  thus  reduce  the  Jeans 
mass  by  a  factor  of  ten.  Because  of  the  cooling  efficiency  of  the 
interstellar  medium  at  low  densities,  the  cooling  time  behind  a  shock 
in  an  KI  region  is  typically  two  orders  of  magnitude  less  than  the  dynamic 
time  scale  (Field  et  al.   1968;  Aanestad  1973),  so  that  the  shock  propagates 
isothermally.  The  jump  in  density  across  an  isothermal  shock  front  is 
approximately  (Kaplan  1966) 

n,      v   , 

tt"   C -tT  ,  (58) 

nl   2  km?1 

where  v  is  the  shock  velocity,  (measured  in  km  s  )  and  may  be  as  large 
as  20  km  s   for  a  strong  shock. 

Various  mechanisms  have  been  proposed  for  producing  and  maintaining 
interstellar  shocks,  and  the  possibility  of  shock- triggered  star  formation 
has  been  examined  under  a  variety  of  physical  conditions.  The  hydro- 
dynamical  models  of  Stone  (1970)  indicate  that  star  formation  may  be 
enhanced  by  shocks  generated  during  collisions  between  interstellar  clouds. 
Indeed,  Loren  (1976)  believes  that  ongoing  star  formation  in  the  NGC  1333 
molecular  cloud  is  the  result  of  such  a  cloud-cloud  collision.  Roberts 
(1969),  Shu  et  al.    (1972),  and  Biermann  et  al.    (1972)  have  suggested  that 
shock  waves  associated  with  density  waves  in  spiral  galaxies  may  induce 
the  gravitational  collapse  of  gas  clouds  thus  leading  to  star  formation. 
Giant  H1I  regions  associated  with  young  early-type  stars  often  line  up 
'like  beads  on  a  string1  along  the  spiral  arms  of  our  Galaxy,  and  there 
is  recent  evidence  for  star  formation  by  density  wave  shocks  in  M33  as 
well  (Dubout-Crillon  1977). 
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The  shock  front  associated  with  the  advancing  ionization  front  of 
an  HEX  region  may  trigger  the  collapse  of  stellar-mass  condensations 
(Dyson  1968) .  Large  OB  associations  may  be  caused  by  a  sequential 
burst  of  HI I  regions  in  a  dense  cloud  (Elmegreen  and  Lada  1977) ,  or 
perhaps  by  a  supernova  cascade  process  (Ogelman  and  Maran  1976) . 
Observations  of  the  Origem  Loop  supernova  remnant  (Berkhvijsen  1974) 
and  the  expansion  of  the  Gum  Nebula  (Schwartz  1977)  suggest  that  the 
strong  shock  from  a  supernova  explosion  may  trigger  star  formation. 
Cameron  and  Truran  (1977)  explain  various  isotoplc  anomalies  and  traces 
of  extinct  ratioactivities  in  solar  system  material  as  being  the  result  of 
a  nearby  Type  II  supernova  that  triggered  the  collapse  of  a  cloud  which 
led  eventually  to  the  formation  of  the  solar  system.  The  detailed  too- 
dimensional  numerical  hydrodynamic  calculations  of  Woodward  (1976) 
demonstrate  the  validity  of  the  shock- induced  mechanism  of  star  formation, 
particularly  when  the  effects  of  self -gravitation,  thermal  instabilities, 
and  dynamical  instabilities  of  the  Kelvin-Helmholtz  and  Rayleigh-Taylor 
type  {cf.   Chandrasekar  1961)  which  are  triggered  by  the  shock,  are  taken 
into  account. 

Thermal  Instabilities 
Thermal  instabilities  in  the  interstellar  medium  can  result  in  the 
fetation  of  non-gravitational  condensations  of  higher  density  and  lower 

temperature  than  are  found  in  the  surrounding  medium  (Field  1965) . 

2 
Basically,  this  is  because  cooling  rates  at  low  densities  vary  as  n 

while  heating  rates  vary  only  as  n.  Following  a  thermal  instability,  as 

the  density  increases  and  the  temperature  drops  (pressure  equilibrium 

obtaining  for  short-wavelength  perturbations) ,  the  critical  Jeans  mass 

given  by  Eq.  (57)  decreases  rapidly.  Theoretical  studies  by  Hunter 


(1966,  1969)  and  Stein  and  McCray  (1972)  have  shown  that  self -gravitating, 
primary,  stellar-mass  condensations  can  form  out  of  the  medium  directly, 
without  the  occurrence  of  fragmentation,  by  the  two-step  process  of 
thermal  instability  at  pressure  equilibrium  followed  by  gravitational 
collapse.  Observationally,  isolated  primary  condensations  having  stellar 
masses  are  known  to  exist  (Aveni  and  Hunter  1967,  1969,  1972;  Herbig  1970, 
1976).  Replacing  the  usual  assumption  of  isothermal  compression  with  the 
condition  of  energy  balance,  Kegel  and  Traving  (1976)  have  generalized  the 

Jeans  criterion  for  gravitational  instability,  and  they  find  that  the 

d£nP  ,,, 
minimum  unstable  mass  is  reduced  by  a  factor  (--^ )   <1,  with  P  and  p 

being  the  pressure  and  density  at  energy  equilibrium. 

Thermal-chemical  instabilities  may  also  lower  the  Jeans  mass.  The 
formation  of  hydrogen  molecules  on  grain  surfaces  in  interstellar  clouds 
may  result  in  pressure  instabilities  leading  to  the  formation  of  protostars 
(Schatzman  1958;  Reddish  1975).  Because  the  cooling  efficiency  is  greater 
for  CO  than  for  CII,  the  conversion  of  CII  to  CO  during  the  evolution  of 
dense  interstellar  clouds  (of.  Herbst  and  Klemperer  1973;  Allen  and 
Robinson  1977)  may  lead  to  instabilities  (Oppenheimer  and  Dalgarno  1975; 
Glassgold  and  hanger  1976).  Generalizing  Field's  (1965)  work  to  include 
chemical  effects,  Glassgold  and  ianger  find  unstable  masses  of  stellar  order. 
Oppenheimer  (1977)  has  demonstrated  that  the  interstellar  gas  may  be 
unstable  to  the  isentropic  growth  of  linear  perturbations  in  dense, 
optically-thick  regions  where  the  molecular  transitions  governing  the 
cooling  of  the  gas  are  thermalized,  and  where  strong  heat  sources  are 
present.  Such  instabilities  may  also  lead  to  the  formation  of  protostars. 

The  criterion  for  thermal  instability  becomes  modified  in  the 
presence  of  magnetic  fields  (Field  1965) .  Just  as  in  the  case  of 
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shock- induced  density  growth  in  a  magnetized  plasma  (of .  Kaplan  1966) , 
magnetic  pressures  inhibit  compression  of  the  gas  in  directions 
perpendicular  to  the  field  lines.  Even  so,  Mufson  (1975)  has  shown  for 
a  wide  variety  of  physical  conditions,  that  the  post-shocked  gas  is 
likely  to  become  thermally  unstable  and  that  condensation  modes  can 
grow  across  magnetic  field  lines.  High  resolution  radio  observations  of 
the  supernova  remnant  IC  443  by  Buin  and  van  der  Laan  (1975)  give  evidence 
for  condensation  perpendicular  to  field  lines. 

Physical  Conditions  in  Dark  Clouds 

Young  stars  (e.g.  T  Tauri  stars,  Herbig  Ae  and  Be  stars,  and  Herbig- 
Haro  objects)  and  pre-stellar  objects  (e.g.  IR  and  maser  sources)  are 
frequently  associated  with  dense  molecular  clouds  (cf.   Strom  et  al .  1975) . 
The  apparent  location  of  newly  formed  stars  and  HII  regions  on  the  outsides 
of  dense,  massive  clouds  and  not  at  their  centers  (Zuckerman  and  Palmer 
1974;  Kutner  et  al.   1976;  Eljnegreen  and  Lada  1977;  Vrba  1977)  suggests  a 
star- formation  scenario  wherein  a  shock-driven  implosion  at  the  boundary  of 
a  cloud  initiates  a  thermal -gravitational  instability,  ultimately  resulting 
in  gravitational ly- bound  condensations. 

Theoretical  studies  by  Solomon  and  Wickramasinghe  (1969)  and  Hollenback 
et  al.   (1971),  and  the  dense  cloud  chemical  models  of  Herbst  and  Klemperer 

(1975)  and  Allen  and  Robinson  (1977) ,  indicate  that  hydrogen  is  pre- 

3   -3 
dominantly  molecular  in  dense  (n  J  10  cm  )  clouds.  Rocket  observations 

by  Carruthers  (1970)  and  Copernicus  satellite  observations  by  Spitzer  et  al. 

(1973)  support  this  conclusion.  Other  major  chemical  constituents  include 

He,  CO,  NH,,  H,  HO,  OH,  HXO,  and  H,0.  A  representative  mean  molecular 

weight  for  dense  cloud  material  would  be  u=2.5.  Dark  clouds  typically 
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3   -3 
have  particle  densities  vpt\,  =  10  cm   and  kinetic  gas  temperatures 

M2 
T=10K  (Heiles  1969;  Penzias  et  al .   1972;  Zuckerman  and  Palmer  1974, 

and  references  cited  therein,  see  also  Appendix  A  of  this  paper  for  a 

detailed  calculation  of  cloud  thermodynamics) . 

Observations  do  not  show  any  interstellar  clouds  rotating  much 

faster  than  the  Galaxy  (cf .   Heiles  1970;  Heiles  and  Katz  1976;  Lada  et  al. 

1974;  Kutner  et  al .   1976;  Bridle  and  Kesteven  1976).  In  the  solar 

-15  -1 
neighborhood,  the  Galaxy  rotates  with  an  angular  velocity  uv,=10    s  . 

Observed  line  widths  of  molecular  transitions  originating  in  dense 
molecular  clouds  are  almost  invariably  too  wide  to  be  explained  by  thermal 
motions,  and  they  frequently  imply  supersonic  velocities.  The  line  widths 
have  commonly  been  attributed  to  microturbulence  (Leung  and  Liszt  1976) 
or  macroturbulence  (Zuckerman  and  Evans  1974),  but  difficulties  with  line 
profile  interpretation  (Snell  and  Loren  1977)  and  energetic  difficulties 
associated  with  supersonic  turbulence  (Dickman  1976,  and  private  communi- 
cation) have  led  to  the  supposition  that  the  line  widths  reflect 
systematic  motions  within  the  clouds,  probably  large-scale  collapse 
(Goldreich  and  Kwan  1974;  Scoville  and  Solomon  1974;  Liszt  et  al.   1974; 
Gerola  and  Sofia  1975;  de  Jong  et  al.   1975;  Snell  and  Loren  1977;  Plambeck 
et  al.   1977;  Fallon  et  al.   1977).  However,  in  favor  of  a  turbulent  origin, 
Arons  and  Max  (1975)  have  suggested  that  the  observed  large  line  widths  may 
be  due  to  the  presence  of  moderate-amplitude  hydromagnetic  waves  in 
molecular  clouds.  Such  waves  may  be  generated  by  the  magnetic  braking 
process . 

Magnetic  field  strengths  in  interstellar  clouds  are  very  uncertain. 
Measuring  Zeeman  splitting  in  the  21  cm  line  of  neutral  hydrogen  and  the 
18  an  OH  line,  Verschuur  (1970)  has  obtained  field  strengths  in  a  number 
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of  diffuse  (n<103  an"3)  clouds.  Clark  and  Johnson  (1974)  have  suggested 
that  the  apparently  anomalous  broadening  of  millimeter-wavelength  SO 
lines  observed  in  Orion  is  caused  by  the  Zeeman  effect  in  6-gauss 
Tuignetic  fields.  However,  Zuckerman  and  Palmer  (1975)  believe  that  the. 
large  line  widths  are  probably  due  to  kinematic  rather  than  magnetic 
broadening.  Beichman  and  Chaisson  (1974)  find  evidence  from  infrared 
polarization  measurements  and  OH  Zeeman  patterns  for  milligauss  fields  in 
the  Orion  infrared  nebula.  Rickard  et  al.   (1975)  and  Lo  et  ai.  (1975) 
have  derived  milligauss  field  strengths  for  a  number  of  OH  maser  sources . 
However,  well-known  observational  and  theoretical  problems  in  interpreting 
Cti  spectra  in  terms  of  Zeeman  patterns  (Zuckerman  and  Palmer  1975;  Heiles 

1976)  make  these  results  tentative.  Magnetic  field  strengths  obtained  by 

3-3  6-3 

Yerschuur  (n<10  cm  ),  Beichman  and  Chaisson  (n=10  cm  ),  and  Lo  et  al. 

8   -3 
(n=10  cm  ),  are  plotted  in  Figure  1  as  a  function  of  inferred  particle 

density  in  the  magnetic  region.  We  employ  Zuckerman  and  Palmer's  estimate 

of  the  gas  density  in  the  Reichman-Chaisson  source,  and  for  the  density 

in  the  neighborhood  of  the  source  discussed  by  Lo  et  al.,   we  take 

u=10  cm   as  suggested  by  Mouschovias  (1976b) . 

2   -3 
At  low  densities  (n<10  cm  ) ,  the  magnetic  field  strength  reflects 

the  large-scale  galactic  field  B„=3uG.  The  constancy  of  the  field 

strength  for  these  low-density  clouds  suggests  that  material  may  stTeam 

preferentially  along  the  field  lines  until  higher  gas  densities  are  reached. 

Anisotropic  gas  flow  along  magnetic  field  lines  increases  the  thermal  gas 

pressure  P  =nkT,  k=1.38xl0"   erg  deg"  being  Boltzmann's  constant,  while 

8  2 

holding  the  magnetic  pressure  P,,=B  /8tt  constant.  Neglecting  inertial 

.  forces,  these  two  pressures  will  come  into  balance  when  the  gas  density 
readies  a  critical  value  given  by 
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ncr  ■  B2/8irkT  ,  (59) 

so  that  for  B=BG=3uG  and  T=10K, 

N   =  260  cm"3  ,  (60) 

which  is  in  good  agreement  with  Figure  1.  Indeed,  the  Parker  (1966) 
instability  (a  magnetic  Rayleigh-Taylor  instability)  may  provide  a 
mechanism  for  preferential  gas  flow  along  magnetic  field  lines  at  low 
densities,  and  the  observational  findings  of  Appenzeller  (1971)  and  Vrba 
(1977)  support  Parker's  predictions. 

Assuming  pressure  equilibrium  maintains  for  n>n  ,  the  magnetic  field 
strength  should  scale  with  the  gas  density  according  to  Eq.  (59): 

B  ■  (8nkT)^  ,  (61) 


so  that 


B  -  nh  (62) 


for  an  isothermal  compression,  in  agreement  with  the  detailed  equilibriun 
models  of  Mouschovias  (1976a,  1976b)  for  self-gravitating  magnetic  clouds. 
This  result  is  to  be  compared  with  Eq.  (1)  which  obtained  for  the  case  of 
isotropic  contraction  and  strict  flux- freezing: 

B  -  n?    ,  (63) 

where  we  have  used 

3 
p  =  until  =  3m/4irR  (64) 

to  relate  the  radius  of  a  spherical  cloud  to  its  particle  density, 
m,  =1.67x10   g  being  the  mass  of  the  hydrogen  atom.  Because  of  the 
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uncertainties  in  determining  B  and  n  for  Figure  1,  it  is  not  possible 
to  determine  precisely  whether  a  slope  of  1/2  or  2/3  best  fits  the 
observations:  neither  is  inconsistent.  However,  as  Mouschovias  (1976b) 
points  out,  a  slope  of  2/3  may  be  incompatible  with  certain  properties 
(e.g.  size,  density,  inferred  magnetic  fields)  of  maser  sources. 
Furthermore,  Scalo  (1977)  has  shown  that  the  heating  of  dense  inter- 
stellar clouds  by  ambipolar  diffusion  imposes  a  constraint  on  cloud  field 
strengths:  B  must  not  increase  faster  than  n  -   so  that  predicted  gas 
temperatures  do  not  exceed  those  observed  in  dense  clouds. 

Initial  Conditions  for  Collapse 
From  an  equation  of  motion  of  the  form  given  by  Eq.  (27) ,  one  can 
derive  [cf.   Cox  and  Giuli  1968)  a  fairly  general  form  of  the  virial 
equation: 

hi  =   2K+3<Y-l>U+M+£2-3PsV  .  (65) 

"     r\2l 

Here  i  =  £_i  >  where  I  is  the  moment  of  inertia  of  the  fluid  about  the 

dt2 
origin  of  coordinates,  K,  U,  M,  and  a   are,  respectively,  the  total  kinetic 

energy  of  mass  motion,  the  thermal  energy,  magnetic  energy,  and  gravi- 
tational energy  within  the  volume  V,  Ps  is  the  hydrostatic  pressure  on 
the  surface  defined  by  V,  and  y   is  the  ratio  of  specific  heats  (y=7/5  for 
a  low- temperature  gas  of  diatomic  molecules) .  From  what  has  been  said 
earlier  regarding  rotation  and  turbulence,  we  may  safely  ignore  the  mass- 
motion  kinetic  energy  term.  Also,  although  strong  surface  pressures  may 
arise  from  passing  shock  waves,  it  is  primarily  the  thermal  instability 
of  the  gas  that  drives  the  condensation  of  material  to  the  higher  densities 
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required  for  eventual  gravitational  collapse.  Accordingly,  we  neglect 

the  P  term  as  well, 
s 

The  condition  for  collapse  is  I<0  .  Eq.  (65)  then  becomes 

|fi|>3<Y-l>U+M  .  (66) 

2 
For  a  uniform  spherical  mass  distribution,  8  "  - '»  -=—  .  Dividing 

Eq.  (66)  by  the  volume  of  the  (spherical)  cloud  V,  noting  that  for  a 

non-relativistic  gas  the  pressure  is  two-thirds  the  energy  density,  and 

assuming  mechanical  equilibrium  between  the  thermal  pressure  (Pp=nkT) 

2 

and  magnetic  pressure  (P.,=B  /8tt)  ,  the  condition  for  gravitational 

collapse,  Eq.  (66),  becomes 

n  >  3.7Sxl05(m/m  )"2  on-3  ,  (67) 

o  w 

or,  equivalently, 

R  <  6.73xl016(m/m_)  cm  ,  (68) 

where  we  have  used  Eq.  (64)  to  eliminate  n  in  favor  of  R,  and  the 
subscripts  here  denote  initial  (i.e.  critical)  values  for  collapse.  The 
initial  magnetic  field  strength  at  thp  cloud  surface  is  found  from 
Eqs.  (59) -(61)  to  be 

B0  =  BG(no/260)!'2  .  (69) 

'Galactic  rotation  (u>r=10    s  )  sets  a  lower  limit  to  the  angular 
velocity  of  a  contracting  cloud,  and  an  upper  limit  is  imposed  by  con- 
servation of  angular  momentum,  provided  there  are  no  external  torques 
acting  on  the  cloud.  Since  the  evolution  of  a  condensation  up  to  the 
time  it  becomes  gravitationally-bound  is  highly  uncertain,  we  do  not 


attempt  to  calculate  the  efficiency  of  magnetic  bTaking  during  this 
stage.  Therefore,  the  angular  velocity  of  a  marginally-unstable  cloud 
cannot  be  determined  a  priori.     It  is  possible  that  a  condensation  may 
derive  its  rotation  from  (subsonic)  turbulence  which  may  be  generated  by 
the  dynamical  instabilities,  particularly  the  Kelvin-Helmholtz  modes 
(of.  Woodward  1976) ,  following  the  passage  of  a  shock.  Because  of  the 
strong  dissipation  of  supersonic  turbulence  (Heisenburg  1948) ,  turbulent 
velocities  must  not  exceed  the  sound  speed 

cQ  -  (YkT/um^  (70) 

-1  7 

which  is  about  0.3  km  s   for  T=10K,  u=2.S,  and  Y=  j  .  If  the  correlation 

length  of  the  turbulence  is  of  the  order  of  the  cloud's  diameter,  then 

u  ■  c  /2R  .  (71) 

o   o  '  ■" 

A  turbulent  origin  for  the  angular  momentum  of  protostars  has  the 
attractive  feature  of  explaining  (1)  the  random  orientation  of  rotational 
axes  of  early-type  stars  (Huang  and  Struve  19S4)  and  field  Ap  stars 
(Abt  et  al,   1972),  (2)  the  lack  of  a  dependence  of  inclinations  in  visual 
binaiy  systems  on  galactic  latitude  (Finsen  1933) ,  and  (3)  the  lack  of 
evidence  (Huang  and  Wade  1966)  for  preferred  galactic  distribution  of 
orientations  of  orbital  planes  of  eclipsing  binaries.  If  stars  and  stellar 
systems  acquired  their  angular  momenta  directly  from  galactic  rotation, 
angular  momentum  vectors  would  generally  be  aligned  perpendicular  to  the 
galactic  plane. 

Initial  values  of  n  ,  R  ,  B  ,  and  u  appear  in  Table  1  for  a  range 
of  protostellar  masses  from  1  m  up  to  40  m   .  The  non-integer  masses 
correspond  to  main-sequence  spectral  types  for  which  main-sequence 
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Table  1.  Initial  values  of  particle  density  n  (cm  ),  cloud  radius 

R  (cm) ,  surface  magnetic  field  strength  B  /Br ,  and  angular 

velocity  u  (s  )  for  various  cloud  masses  m/m  marginally 

unstable  to  gravitational  collapse.  Numbers  in  parentheses 
are  decimal  exponents. 


m/m  Spectral  n  (cm"3)  R  (cm)  B  /B„  u  (s-1) 

e  r,  o  o  ^  '  oG  o 
type 

1  G2  3.75  (5)  6.73  (16)  38.0  2.23  (-13) 

1.7  F0  1.30  (5)  1.14  (17)  22.4  1.31  (-13) 

2.1  AS  8.50  (4)  1.41  (17)  18.1  1.06  (-13) 

3.24  A0  3.57  (4)   ■  2.18  (17)  11.7  6.88  (-14) 

6.5  B5  8.88  (3)  4.37  (17)  5.8  3.43  (-14) 

10  B3  3. 75  (3)  6.73  (17)  3.8  2.23  (-14) 

17.8  B0  1.18  (3)  1.20  (18)  2.1  1.24  (-14) 

40  05  2.34  (2)  2.69  (18)    1  5.57  (-15) 
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rotation  data  have  been  accumulated.  Spectral  types  later  than  early  G 

will  not  be  considered  since  these  stars  are  expected  to  lose  most  of 

their  primordial  angular  momentum  during  pre-main-sequence  contraction 

and  ma in- sequence  nuclear  burning  (cf.   Section  I).  Interestingly  enough, 

the  values  of  u    are  roughly  the  same  as  those  one  would  calculate  assuming 

conservation  of  angular  momentum  for  a  condensation  having  initial 

densities  of  the  order  of  a  few  particles  per  cubic  centimeter,  and 

initially  rotating  with  the  Galaxy.  Clearly  then,  the  adopted  values  for 

u  are  probably  overestimated.  By  adopting  a  possibly  exaggerated  oiQ,  we 

require  the  magnetic  braking  to  be  correspondingly,  more   efficient  in 

decelerating  a  protostar.  (It  will  turn  out  that  the  calculations  are 

quite  insensitive  to  a  wide  range  of  values  of  o)q.)  The  values  of  nQ 

and  01  for  the  1  m  and  2.1  m  clouds  are  in  agreement  (fortuitously) 
o         ©        © 

with  the  initial  conditions  assumed  by  Black  and  Bodenheimer  (1976) 

in  their  calculations  of  rotating  protostars.  For  all  mases,  the  initial 

3  2 
ratio  of  gravitational  to  centrifugal  forces  at  the  equator,  Fg/Fc=Gm/R0"0. 

is  about  ten.  Thus,  centrifugal  forces  are  not  sufficient  to  stabilize 

the  initial  configurations,  and  the  assumption  that  2K«fi  appears 

reasonable.* 


Goldreich  and  Lynden-Bell  (1965)  and  Toomre  (1964)  have  shown  by 
a  generalization  of  the  Jeans  stability  criterion  to  include  the  sta- 
bilizing effect  of  rotation,  that  whereas  in  the  classical  Jeans  case 
with  pressure  effects  stabilizing  short -wavelength  perturbations,  long 
waves  are  stabilized  by  rotation. 


SECTION  IV 
MAGNETIC  BRAKING  OF  COLLAPSING  PROTOSTARS 

We  now  proceed  to  calculate  the  rotational  deceleration  of  a 
magnetically  braked,  contracting  protostar.  The  value  of  the  rotational 
deceleration  is  given  by 

&>  _   2R     t  ,,7, 

dt     R     KMR *  . 

With  the  initial  conditions  given  in  Table  1,  u>=u)(R)  is  obtained  by 
numerical  integration  of  Eq.  (47)  together  with  R,  t,  and  B  (R)  given  by 
Eqs.  (41),  (18),  and  (44),  respectively.  A  fourth-order  Runga-Kutta 
integration  scheme  with  variable  step-size  was  employed.  The  accuracy 
of  the  numerical  code  was  tested  by  setting  t=0  in  Eq.  (47)  and  checking 
conservation  of  angular  momentum  (J=0) .  The  value  of  the  cloud  radius 
at  each  step  in  the  integration  was  computed  by  a  subroutine  which  solved 
Kepler's  equation  for  a  collision  orbit  (i.e.,  a  degenerate  ellipse  with 
eccentricity  e=l)  by  a  standard  iterative  procedure.  This  was  found  to  be 
easier  than  integrating  Eq.  (41)  and  solving  the  resulting  transcendental 
equation  for  R  at  each  step. 

In  oider  to  keep  the  problem  tractable,  a  number  of  simplifying, 
although  reasonable,  assumptions  have  been  made.  To  simplify  the  geometry, 
the  collapse  is  assumed  to  be  isotropic  and  homologous,  which  implies 
further,  that  the  surface  poloidal  magnetic  field  increases  according  to 

Bs(R)  =  BQ(Ro/R)2  ,  (72) 
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with  B  taken  from  Table  1.  Furthermore,  it  is  assumed  that  the  magnetic 
stresses  within  a  cloud  will  constrain  the  cloud  to  rotate  uniformly 
as  a  rigid  body.  It  may  be  argued  that  the  assumption  of  homologous 
collapse  may  be  somewhat  artificial  in  view  of  the  numerical  collapse 
models  of  Larson  (1969,  1972b)  and  Hunter  (1969).  However,  there  is  little 
or  no  observational  evidence  for  'Larson-type'  dynamical  evolution 
(Cohen  and  Kuhi  1976) ,  and  Disney  T1976)  believes  the  impressed  boundary 
conditions  in  the  Larson  approximation  are  probably  not  realistic. 
A  nonhomologous  collapse  would  have  the  effect  of  lowering  the  moment 
of  inertia  of  a  contracting  cloud,  as  well  as  increasing  the  gravitational 
(binding)  energy  somewhat. 

The  initial  value  of  the  azimnthal  magnetic  field  B, (RQ)  can  not  be 
determined  a  priori.  We  adopt  B  (Rq)=0,  with  the  understanding  that  the 
braking  efficiency  will  be  somewhat  underestimated   during  the  initial 
collapse  phase  since  x  =0. 

Because  the  magnetic  field  of  a  cloud  remains  frozen- in  during  the 
collapse,  the  free-fall  time  given  by  Eq.  (56)  which  is  valid  only  for 
a  pressure-free  collapse,  will  underestimate  the  true  collapse  time. 
From  the  equilibrium  virial  theorem  [1=0  in  Eq.  (65)]  ,  one  can  define  an 
effective  Newtonian  gravitational  constant 

,  „    |n|-2K-M-3<yl>U:)   _  (73) 

|a| 

Since  the  collapse  proceeds  isothermally  at  the  relatively  low  densities 
(n<108  cm"3)  considered  here  (Gaustad  1963;  Gould  1964;  Hayashi  and  Nakano 
1965;  Hattori  et  al.   1969;  see  also  Appendix  A  of  this  paper),  the  term 
3<y-1>U  representing  thermal  gas  pressure  is  not  expected  to  contribute 
significantly  to  Eq.  (73).  Furthermore,  provided  a  cloud  loses  angular 
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momentum  during  collapse,  the  term  2K  can  be  neglected  (see  discussion 
in  Section  III).  Thus,  the  collapse  of  a  magnetic  protostar  proceeds 
on  a  magnetically-diluted  time  scale  given  by 

t  =  (JUL   Vs  (74) 

where  B'  »  G(]fi|-H)/|n|  .  Since  |Q J  -  2MQ  (see  Eq.  (66)  and  discussion 
leading  to  Eq.  (67)),  and  both  |«|  and  M  grow  like  ~R      for  an  isotropic 
collapse  with  flux  conservation,  G'=  G/2  . 

The  magnetic  braking  mechanism  is  operative  only  so  long  as  a  cloud 
remains  magnetically  linked  to  its  surroundings.  Mestel  (1966)  and 
Mestel  and  Strittmatter  (1967)  have  argued  that,  as  a  cloud  contracts, 
the  almost  oppositely  directed  field  lines  at  the  equatorial  plane  give 
rise  to  strong  'pinching'  forces  that  dissipate  flux  and  reconnect  field 
lines,  so  that  the  magnetic  field  of  a  cloud  is  effectively  detached  from 
that  of  the  background.*  However,  as  Mestel  and  Strittmatter  themselves 
point  out,  the  time  scale  for  this  process  is  so  long  that  this  process  may 
not  be  efficient.  Furthermore,  the  equilibrium  models  of  Mouschovias  (1976b) 
do  not  show  any  tendency  for  equatorial  pinching. 

The  expulsion  of  a  cloud's  magnetic  field  by  ambipolar  diffusion 
provides  a  more  efficient  mechanism  for  detaching  the  cloud's  field  from 


Kulsrud's  (1971)  mechanism  for  the  rotational  deceleration  of  a 
rotating  dipole  applies  in  this  case.  Kulsrud  finds  t~FPBsB0.  However, 
because  the  magnetic  torques  are  then  proprotional  to  the  magnetic  field 
in  the  surrounding  medium  instead  of  the  azimuthal  field,  the  braking 
efficiency  is  reduced  considerably.  Kulsrud's  formulae,  derived  for  a 
harmonically  time-varying  dipole,  may  not  be  strictly  applicable  to  a 
contracting  protostar,  the  radius  of  which  decreases  secularly  with 
time. 
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that  of  its  surroundings  (Mestel  and  Spitzer  1956;  Nakano  and  Tadomaru 

1972) .  The  time  scale  for  ambipolar  diffusion  is  (Nakano  and  Tademaru 

1972) 

?  -2 

t„  =  87iRmii<av>nn  B 
D  "2  e 

=  8.26xl021   (ne/n)      ,  (75) 

where  we  have  used  Eqs.  (64),  (67) -(69),  and  (72)  to  write  the  second 

-9   3  -1 
equality,  taking  <av>=2xl0   cm  s   (Osterbrock  1961).  The  electron 

density  n  is  calculated  from  Eq.  (B7)  which  appears  in  Appendix  B 

(see  Appendix  B  for  a  discussion  on  the  fractional  ionization  in  dense 

magnetic  clouds) .  Because  the  ionozation  rate  in  dense  clouds  is  somewhat 

uncertain  (cf.   Appendix  B) ,  we  consider  the  two  limiting  cases:  weak 

40 

ionization  by  radioactive   K  nuclei  only,  at  a  rate  given  by  £„ 

(Eq.  (B9)),  and  ionization  by  both   K  and  (magnetically  screened) 
cosmic  rays  at  a  rate  determined  essentially  by  C™  (Eq.  (B19)),  bearing 
in  mind  that  the   K  rate  is  probably  the  more  realistic  of  the  two 
(cf .  Appendix  B) .  As  pointed  out  by  Mouschovias  (1977)  ,  if  there  is 
tension  in  the  field  lines,  Eq.  (75)  overestimates  the  diffusion  time 
scale.  However,  the  field  inside  the  cloud  is  assumed  uniform  so  that 
tj,  is  probably  not  much  less  than  that  given  here  by  Eq.  (75). 

The  magnetic  field  becomes  essentially  uncoupled  when  t„  -   t 
(Mestel  and  Spitzer  1956;  Nakano  and  Tademaru  1972).  Afterwards, 
a  cloud  contracts  conserving  angular  momentum.  Thus,  the  angular 
momentum  of  a  protostar  is  established  at  the  uncoupling  epoch.  The 
integrations  were  therefore  terminated  when  this  condition  was  met.  For 
a  uniform,  spherical,  rigidly-rotating  protostar,  the  angular  momentum 
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at  the  uncoupling  epoch  is  simply 

Ju  -  0.4  mR^u  .  (76) 

Results  of  the  calculations  appear  in  Table  2. 

In  all  cases  considered  here,  the  magnetic  torques  rotationally 
decelerate  the  clouds,  constraining  them  to  co-rotate  with  their 
surroundings  at  an  angular  velocity  u>  =w  =10"1S  s"1.  Physically  the 

U   b  *  '  * 

Alfven  speed  just  outside  a  cloud  is  always  greater  than  the  collapse 
velocity  at  the  cloud  surface,  so  that  the  torques  are  able  to  transmit 
angular  momentum  from  a  collapsing  cloud  on  a  time  scale  which  is  less 
than  the  (magnetically-diluted)  free-fall  time.  Thus,  the  supposition 
that  magnetic  braking  is  efficient   all    the  Way   down   to   the  breakdown  of 
flux-freezing   (e.g.  Hoyle  1960;  Mouschovias  1977)  appears  to  be  vindicated. 
The  angular  momentum  of  an  initial  condensation  is  reduced  some  three 
to  four  orders  of  magnitude  for  the  case  of  ionization  by  40K  and  cosmic 
rays,  respectively.  The  angular  momentum  of  the  40K  clouds  at  the 
uncoupling  epoch  is  an  order  of  magnitude  greater  than  that  of  the  cosmic- 
ray  ionized  clouds  because  uncoupling  occurs  earlier  in  the  collapse 
sequence  (nu=106  cm"3)  for  the  weaker  40K  ionization  rate.  Clouds  ionized 
primarily  by  cosmic  rays  uncouple  from  their  surroundings  later  (n  =107  cm"3) 
since  the  time  scale  for  ambipolar  diffusion  is  longer  for  these  clouds. 

The  surface  magnetic  field  strengths  at  the  uncoupling  epoch  fall 
nicely  within  the  range  of  those  observed  in  dense  clouds  (see  discussion 
in  Section  III  leading  up  to  Figure  1).  The  ratio  |B  /Bj  at  uncoupling 
was  found  to  be  near  unity  for  all  clouds.  This  is  consistent  with  the 
findings  of  Gillis  et  al.    (1974)  who  predict  a  similar  ratio  in  their  time- 
independent  magnetic  braking  model.  If  |B  /Bj  become  much  larger  than 
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assuming  rigid  rotation,  and  an  upper  bound  is  estimated  by  adopting 

Bodenheimer  and  ©striker's  (1970)  differential  rotation  law  for  the 

early- type  stars  (see  Appendix  C) .  For  rapidly  rotating  stars,  gravity 

darkening  may  lead  to  an  underestimate  of  the  observed  equatorial 

velocity  by  as  much  as  40  percent  (Hardorp  and  Strittmatter  1968) ,  so 

that  the  transition  of  j  between  the  early-type  main-sequence  stars  and 

the  eclipsing  binaries  is  relatively  smooth. 

The  decline  in  j  for  stars  later  than  spectral  type  early-F 

(m/m  <  2)  is  thought  to  be  the  result  of  angular  momentum  transfer  during 
& 

pre-main-sequence  contraction  (Schatzman  1962;  Mestel  1968;  Schwartz 
and  Schubert  1969;  Okamoto  1969,  1970)  and  main- sequence  nuclear  burning 
(Dicke  1964;  Brandt  1966;  Nfodisette  1967;  Weber  and  Davis  1967),  or 
perhaps  an  indication  of  the  presence  of  planetary  systems  (Hoyle  1960; 
McNally  1965;  Huang  1973,  and  references  cited  therein). '  Tarafdar  and 
Vardya  (1971)  account  for  this  discrepancy  in  j  by  presuming  a  rapidly 
rotating  interior  for  the  later- type  stars  and/or  a  slowly  rotating 
interior  for  the  early- type  stars. 

The  important  point  illustrated  by  Figure  2   is  that  if  angular 
momentum  is  conserved  during  star  formation,  then  the  angular  momentum 
of  a  protostellar  condensation,  being  almost  two  orders  of  magnitude 
greater  than  that  of  the  widest  separated  visual  binaries,  will  be  incon- 
sistent with  the  observed  angular  momenta  of  stellar  systems.  The 
two-dimensional  numerical  hydrodynamic  models  of  Larson  (1972a)  and  Black 
and  Bodenheimer  (1976)  for  rotating,  collapsing  protostars  are  thus  highly 
suspect  and  probably  not  physically  realistic  since  they  assume  strict 
angular  momantum  conservation  throughout  the  collapse.  Even  if  the 
toroidal  figures  predicted  by  their  models  (for  which  there  is  no 
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observational  evidence)  are  unstable  to  nonaxi symmetric  breakup,  as 
suggested  by  Wong's  (1974)  stability  analysis  of  equilibrium  toroids, 
the  angular  momentum  of  the  system  remains  unchanged,  and  one  is  hard 
pressed  to  find  a  mechanism  to  dispose  of  angular  momentum. 

On  the  other  hand,  the  calculations  presented  here  for  magnetic 
braking  during  star  formation  are  consistent  with  the  observations 
presented  in  Figure  2.  Single  stars  are  rare  (Blaauw  1961;  Heintz  1969; 

Abt  and  Levy  1976) .  It  is  therefore  not  surprising  that  the  most  likely 
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ionization  rate  in  dense  clouds,  namely,  ionization  by   K  radioactive 

nuclei  (cf .  Appendix  B) ,  predicts  angular  momenta  corresponding  to  that 
observed  for  close  binaries,  while  the  much  less  likely  situation  of 
ionization  by  cosmic  rays  accounts  nicely  for  the  angular  momenta  of 
(rare)  single  stars.  The  minimum  angular  momenta  of  single  main-sequence 
stars  is  in  excellent  agreement  with  the  minimum  angular  momentum  of  a 
contracting  protostar,  jfR  .  Single  main-sequence  stars  with  rapidly _ 
rotating  cores  and/or  large  equatorial  velocities  apparently  form  in 
regions  of  lower  cosmic  ray  flux.  These  are  the  stars  in  Figure  2  having 

JCR<j<V 

Apparently,  clouds  having  a  specific  angular  momenta  JK  become 

unstable  ana  fragment  into  a  multiple  (e.g.  binary)  star  system.  After 

the  magnetic  field  is  expelled  from  a  contracting  cloud,  the  cloud 

continues  to  contract  conserving  angular  momentum  since  no  external 

(magnetic)  torques  act  on  the  cloud.  Eventually,  as  the  rotational 

kinetic  energy  of  the  cloud  increases,  the  ratio  of  centrifugal  forces  to 

gravity  exceeds  a  critical  value  determined  essentially  by  the  distribution 

of  mass  and  angular  momentum  within  the  configuration  (see  Ostriker  1970 

for  references) ,  and  the  cloud  breaks  up  into  two  or  more  condensations 
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unity,  there  would  be  the  danger  that  hydro-magnetic  instabilities  of 
the  twisted  field  might  interfere  with  the  assumed  poloidal  magnetic 

topography. 

3  2 

The  ratio  F  /F  =  Gm/R  id  indicates  the  important  result  that 
g  c      u  u 

centrifugal  forces  are  kept  well  below  gravity   throughout  the  collapse 
sequence. 


SECTION  V 
DISCUSSION  OF  RESULTS:  COMPARISON  WITH  OBSERVATIONS 

The  brightest  flashes  in  the  world  of  thought  are 
incomplete  until  they  have  been  proven  to  have  their 
counterparts  in  the  world  of  fact. 

-  John  Tyndall 

(British  Physicist  1820-1893) 

The  hypothesis  of  magnetic  braking  during  star  formation  offers  a 
plausible  explanation  for  the  observational-  fact  that  interstellar  clouds, 
in  general,  do  not  rotate  much  faster  than  the  Galaxy  (Heiles  1970; 
Heiles  and  Katz  1976;  Bridle  and  Kesteven  1976;  Kutner  et  al .   1976; 
Loren  1977  and  private  communication;  Lada  et  al.   1974).  Further  obser- 
vational evidence  for  angular  momentum  transfer  during  star  formation  is 
afforded  by  a  consideration  of  the  angular  momenta  of  binary  systems 
and  single  main-sequence  stars,  stellar  rotation  on  the  main  sequence, 
and  the  angular  momentum  of  the  protosun. 

Specific  Angular  Momenta  of  Single  and  Binary  Stars 
The  specific  angular  momenta  j  (angular  momentum  per  unit  mass)  for 
binary  systems  (visual,  spectroscopic,  and  eclipsing),  single  main-sequence 
stars ,  and  the  solar  system  (SS) ,  is  illustrated  in  Figure  2 .  In 
calculating  j  for  the  main-sequence  stars,  we  have  assumed  that  the  stars 
have  a  mass  distribution  given  by  Eddington's  standard  model  (polytropic 
index  n  =  3) .  Taking  the  observed  mean  equatorial  rotational  velocities 
for  main-sequence  stars  (Allen,  1973) ,  a  lower  bound  for  j  is  obtained 
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with  most  of  the  angular  momentum  of  the  original  cloud  going  into 
orbital  motion  of  the  fission  fragments.  The  details  of  the  'fission' 
■process  have  yet  to  be  worked  out. 

That  a  cloud  having  a  specific  angular  momentum  j.,  is  expected 
to  fission  can  be  seen  most  easily  in  Figure  3.  Here  we  plot  the 
specific  angular  momenta  for  a  number  of  eclipsing  binary  systems  for 
which  absolute  orbital  dimensions  have  been  determined.  Most  of  the 
data  are  taken  from  Kopal  (1959) .  Roxburgh  (1966)  and  Bodenheimer  and 
Ostriker  (1970)  have  calculated  the  threshold  (i.e.  minimum)  angular 
momentum  necessary  for  fission  to  occur;  Roxburgh  (R)  for  the  W  Ursae 
Majoris  systems,  and  Bodenheimer  and  Ostriker  (BO)  for  early- type  close 
binary  systems.  Their  results,  reproduced  in  Figure  3,  are  shown  to  be 
in  good  agreement  with  the  specific  angular  momenta  of  close  binary  systems, 
and  more  importantly  here,  indicate  that  our  primary  condensations, 
having  a  specific  angular  momentum  j„,  will  be  unstable  to  bifurcation 
since  j„  is  greater  than  the  threshold  j  for  all  masses. 

According  to  our  theory  of  magnetic  braking,  j„  is  an  upper  limit 
to  the  specific  angular  momentum  of  a  cloud  of  given  mass.  How  then  are 
the  wide  (i.e.  long-period  spectroscopic  and  visual)  binaries  formed? 
Evidently,  an  independent  mode  of  formation  exists  for  these  systems. 
Indeed,  there  has  been  increasing  evidence  for  two  separate  modes  of 
binary  formation.  Contrary  to  the  earlier  suggestions  of  Kuiper  (19S5) , 
Van  Albada  (1968a)  finds  that  the  division  of  early-type  binaries  into 
close  (spectroscopic)  and  wide  (visual)  pairs  is  probably  real  and  not  due 
to  the  obvious  selection  effects.  Whatever  the  period  distribution,  smooth 
or  bimodal,  Huang  (private  communication  to  H.A.  Abt  and  S.G.  Levy,  1976) 
feels  that  no  single  formation  process  will  produce  binaries  with  such  a 
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wide  range  of  periods  (1  day<P <  10  days).  From  a  statistical  analysis 
of  the  frequency  of  binary  secondary  masses,  Abt  and  Lev)'  (1976)  conclude 
that  there  are  two  types  of  binaries:,  those  with  the  shorter  periods  are 
fission  systems  in  which  a  single  protostar  subdivided  because  of 
excessive  angular  momentum,  whereas  the  longer  periods  represent  pairs 
of  protostars  that  contracted  separately  but  as  a  common  gravitationally- 
bound  system.  This  'neighboring-condensation'  or  'early- capture1  mechanism 
for  the  origin  of  long-period  binaries  is  supported  by  the  numerical 
calculations  of  Van  Albada  (1968b)  and  Amy  and  Weissman  (1973) ,  as  well 
as  by  the  findings  of  the  recent  two-body  tidal  capture  theory  of  Fabian 
et  al.   (1975;  see  also  Press  and  Teukolsky  1977). 

It  is  not  likely  that  a  wide  binary  system  can  evolve  into  a  close 
system  by  disposing  angular  momentum,  thereby  making  unnecessary  a 
separate  formation  mechanism  for  close  binary  systems.  Binary  stellar 
winds  (Mestel  1968;  Siscoe  and  Heinemann  1974)  can  operate  only  in  the 
later-type,  low-mass  stars  which  have  outer  convection  zones.  Angular 
momentum  loss  via  gravitational  radiation  is  efficient  only  for  low-mass 
systems  already  in  near  contact  (Webbink  1976) ,  and  the  disposal  of 
angular  momentum  by  mass-loss  is  not  expected  to  occur  until  the  late 
main- sequence  evolutionary  stages  of  already  close  binary  systems 
(Webbink  1976).  Furthermore,  Webbink  (1977)  believes  that  most  W  Ursae 
Majoris  systems  have  always  existed  in  a  contact  state,  and  that  fission 
is  the  only  obvious  formation  mechanism  satisfying  this  requirement. 

Rotation  of  Main -Sequence  Stars 
The  equatorial  rotational  velocity  of  a  main-sequence  star  having  a 
specific  angular  momentum  j  is 

v*  "  J/kR*  (77) 
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where  tc  is  the  effective  gyration  constant  for  the  star  having  a  main- 
sequence  radius  R*.  Assuming  angular  momentum  is  conserved  after  the 
uncoupling  epoch,  j  is  given  by  j,.  or  j_.  Since  stars  later  than  early-F 
may  lose  large  quantities  of  angular  momentum  during  pre-mairi-sequence 
contraction,  we  consider  only  the  eaily-type  (m/m_  >  2)  stars.  The  mass 
distribution  within  such  stars  is  given  to  a  good  approximation  by  the 
Eddington  standard  model,  characterized  by  a  polytropic  index  n=3.  For 
rigid  rotation  k=0.08,  and  for  the  differentially  rotating  models  of 
Bodenheimer  and  Ostriker  (1970),  k=0.2.8  (see  Appendix  C  for  details  of 
these  calculations).  Using  this  information  together  with  the  angular 
momentum  data  in  Table  2,  rotational  velocities  are  calculated  from  Eq.  (77) 

and  the  results  appear  in  Table  3.  Independent  of  the  assumed  rotation 

40 
law,  v*  for  the   K  ionization  rate,  far  exceeds  the  critical  equatorial 

breakup  velocities  for  all  spectral  types  (Slettebak  1966)  ;  a  fortiori 
these  stars  are  expected  to  fission  into  a  close  pair.  Within  the  frame-  . 
work  of  our  present  theory  of  magnetic  braking,  single  stars  are  believed 
to  acquire  an  amount  of  angular  momentum  determined  by  the  cosmic  ray 
ionization  rate.  Figure  4  illustrates  the  reasonable  agreement  parti- 
cularly for  the  case  of  uniform  rotation  (k=0.08),  of  predicted  rotational 
velocities  with  those  observed  for  single  main-sequence  stars  (Allen 
1973).  Rigid  rotation  may  result  from  the  actions  of  circulation  currents, 
convective  mixing  in  the  core,  or  magnetic  viscosity  (Roxburgh  and 
Strittmatter  1966),  so  that  the  case  k=0.08  may  in  fact  represent  the 
most  plausible  sittiation.  The  anomalously  high  rotational  velocities 
predicted  for  the  lower-mass  stars  on  the  assumption  of  angular  momentum 
conservation  during  contraction  to  the  main  sequence,  is  evidence  for 
further  rotational  braking  during  the  later  pre-main-sequcnce.  evolutionary 
stages  of  these  stars. 


ss 


Table  3.  Predicted  equatorial  rotational  velocities  v*(km  s  )  for 
uniformly  rotating  (k=0.08)  and  differentially  rotating 
(k=0.28)  main-sequence  stars  for  the  two  limiting  ionization 
rates,  CK  and  £_R.  The  last  column,  taken  from  Allen  (1973), 

gives  mean  values  for  observed  stars. 
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k         - 
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G2 
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710 
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1.7 
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44 
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Angular  Momentum  of  the  Protosun 

Hoyle  (I960,  1963)  has  estimated  the  angular  momentum  of  the  early 

51     2-1 
solar  system  J£SS  =  4x10   g  an  s   by  augmenting  the  planets  up  to 

normal  solar  composition.  This  value  is  in  good  agreement  with  our 

51    2   1 
predicted  Ju  ■  3x10   g  cm  s   (see  Table  2)  for  a  one- solar-mass  cloud 

ionized  primarily  by   K. 

The  specific  angular  momentum  of  the  early  solar  system  (ESS)  is 

compared  with  that  of  the  present  solar  system  (SS)  in  Figure  3.  The 

fact  that  JESS>  lying  well  above  the  threshold  j  necessary  for  fission, 

is  comparable  to  the  specific  angular  momenta  of  close  binary  systems, 

suggests  that  both  (close)  binary  and  planetary  systems  may  be  formed 

by  a  similar  process  (cf.  van  den  Heuvel  1966;  Fleck  1977)  involving 

the  rotational  instability  of  a  single  primary  condensation.*  Drobyshevski 

(1974)  has  suggested  a  mechanism  for  close-binary  formation  wherein  the 

convective  outer  layers  of  a  rapidly  rotating  protostar  are  thrown  off 

forming  a  ring  in  the  star's  equatorial  plane  which  becomes  unstable  and 

forms  a  second  component.  This  is  very  similar  to  the  generally  accepted 

Kant-Laplace  nebular  hypothesis  for  solar-system  formation.  The  criteria 

for  determining  whether  the  end-product  of  such  an  instability  will  be 

planetary  of  stellar  would  be  of  interest.  Circumstellar  disks  commonly 


* 
In  the  past  (cf.   Brosche  1962;  McNally  1965),  the  angular  momentum 
presently   m  the  solar  system  (SS)  had  always  been  compared  with  that  of 
single  mam- sequence  stars.  Such  a  comparison  is  probably  fortuitous 
because,  as  can  be  seen  from  Figure  5,  a  condensation  having  a  specific 
angular  momentum  jgs  is  not  expected  to  become  rotationally  unstable.  One 
must  then  look  for  a  mechanism,  other  than  the  popular  nebular  hypothesis 
wherein  a  planetary  system  forms  in  the  equatorial  plane  of  a  rotationally 
unstable  protostar,  to  account  for  the  formation  of  a  planetary  system. 
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associated  with  Be  stars  may  be  the  result  of  a  rotational  instability, 
particularly  since  these  stars  rotate  with  near-breakup  velocities 
(Slettebak  1966). 


SECTION  VI 
CONCLUDING  REMARKS 

Without  rotational  braking  during  star  formation,  single  stars 
would  rotate  with  speeds  close  to  that  of  light.  This  statement  is 
actually  a  reductio  ad  absurdum;   centrifugal  forces  will  halt  the  col- 
lapse perpendicular  to  the  rotation  axis  long  before  such  speeds  are 
attained.  However,  even  if  the  resulting  highly  flattened  system  frag- 
ments with  most  of  the  angular  momentum  going  into  orbital  motion  of  the 
fragments  about  their  center  of  mass,  multiple  star  systems  (e.g.  binary 
stars)  would  have  periods  two  to  three  orders  of  magnitude  longer  than 
the  periods  typically  observed  for  even  the  widest  pairs. 

The  present  investigation  suggests  that  magnetic  torques  acting  on 
a  rotating,  contracting  cloud  which  is  permeated  by  a  frozen- in  magnetic 
field  coupling  the  cloud  to  its  surroundings,  rotationally  decelerate  a 
cloud,  constraining  it  to  co-rotate  with  the  background  medium.  Centrifugal 
forces  are  always  kept  well  below  gravity.  The  angular  momentum  of 
magnetically  braked  clouds  is  consistent  with  the  observed  angular  momenta 
of  close  binary  systems  and  single  early-type  main-sequence  stars.  The 
hypothesis  of  magnetic  braking  offers  support  to  the  fission  theory  for 
the  formation  of  close  binary  systems,  and  is  able  to  account  for  the 
relative  paucity  of  single  stars.  The  calculations  also  suggest  a  common 
mode  of  formation  for  (close)  binary. and  planetary  systems. 

Throughout  this  investigation,  some  simplifying  assumptions  (most 
of  which  are  physically  justifiable)  have  been  made  in  order  to  keep  the 
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problem  tractable.  In  every  instance,  we  have  deliberately  underestimated 
the  efficiency  of  magnetic  braking.  Even  so,  the  braking  is  still  able 
to  constrain  the  clouds  to  co-rotate  with  their  surroundings;  a  fortiori 
a  more  accurate  analysis  will  reach  the  same  conclusion.  The  one  single 
factor  most  important  in   determining   the  angular  momentum  of  a  protostar 
is   the  ionization  rate  in   dense  magnetic  clouds:    the  degree  of  ionization 
controls  the  coupling  of  a  cloud  to  the  galactic  magnetic  field.  For  this 
reason,  the  material  discussed  in  Appendix  B  concerning  the  fractional 
ionization  in  dense  clouds  is  in  need  of  further  study,  particularly  the 
mechanisms  for  excluding  cosmic  rays  from  magnetic  clouds. 

Our  calculations  suggest  that  the  toroidal  configurations  obtained 
in  the  numerical  collapse  models  of  Larson  (1972a)  and  Black  and 
Bodenheimer  (1976) ,  wherein  angular  momentum  is  conserved  throughout  the 
collapse,  may  not  appear  until  the  later  pre-main-sequence  evolution  of  a 
protostar.  Star  formation  is  a  complicated  process  and  the  (simplifying) 
assumption  commonly  made  that  magnetic  fields  (and  therefore  their  role  in 
magnetic  braking)  can  be  ignored  may  not  be  physically  realistic.  Indeed, 
it  may  be  safer  to  ignore  effects  of  rotation  on  the  evolution  of  a 
protostar!  The  work  of  Larson  and  Black  and  Brodenheimer  should  be 
re-evaluated  incorporating  the  possible  effects  that  magnetic  fields  may 
have  on  the  structure  and  evolution  of  a  protostar. 


APPENDIX  A 
HEATING  AND  COOLING  RATES  IN  DENSE  CLOUDS 

Heating 
Cosmic-ray  heating  and  compressional  heat  generated  by  the  collapse 
are  the  dominant  heating  mechanisms  in  dense  molecular  clouds.  Heating 
by  photo-dissociation  of  H~  (Stephens  and  Dalgamo  1973) ,  by  photo- 
electrons  ejected  from  grains  (Spitzer  1948),  by  photoionization  of  the 
gas  (Takayanagi  and  Nishimura  1960) ,  and  by  chemical  reactions  (Dalgarno 
and  Oppenheimer  1974)  is  unimportant  in  dense  clouds  because  the  ultra- 
violet photons  are  mostly  screened  out.  Heat  of  formation  released  by 
newly  formed  H7  molecular  (Spitzer  and  Cochran  1973)  is  unimportant 
because  of  the  low  neutral  hydrogen  abundance  in  dense  clouds.  Scalo 
(1977)  has  suggested  that  the  action  of  ambipolar  diffusion  may  generate 
an  appreciable  amount  of  heat  if  the  magnetic  field  in  a  contracting  cloud 
grows  like  B~n  where  k>-r  .  This  mechanism  is  probably  unimportant  for 

our  clouds  which  are  characterized  by  k  =  j    during  the  initial  compression 

2 
stages  and  k  =  ^  for  gravitational  collapse.  Heating  by  the  dissipation 

of  (supersonic)  turbulence  is  ignored. 

The  cosmic-ray  heating  rate  per  unit  volume  is 

The  cosmic-ray  ionization  rate  per  H  molecule  £_.  is  computed  in  Appendix  B 
and  is  given  by  Eq.  (B19) .  Glassgold  and  Langer  (1973a)  give  the  mean 
energy  gain  per  ionization  <E.  >=17  eV.  The  heating  rate  by  freely 
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propagating  cosmic  rays  (f=l  in  Eq.  (B19))  is  then 


'CR 


1  2 

2.72xl0~"  n  exp[- 1.58x10  (m/m  )  n  ]  erg  s"1  cm   .    (A2) 


A  measure  of  the  heat  generated  by  the  collapse  r  is  the  ratio  of  the 
thermal  energy  density  of  a  cloud  at  temperature  T,  j  nkT  (k=l. 38x10 
erg  deg    is  Boltzmann's  constant),  to  its  free-fall  time  tf  defined 
by  F.q,  (56).  Thus 

r  «  2xl0"31ri"2"  T  erg  s"1  cm"3  .  (A3) 


Molecular  Cooling 

Inelastic  gas-grain  collisions,  and  rotational  transitions  among 
the  more  abundant  molecular  species,  H-,  CO,  and  IID,  will  cool  the  gas. 
Other  perspective  molecular  coolants  such  as  H,CO  (Thaddeus  1972) ,  HC1 
(Dalgamo  et  al .   1974) ,  and  CS  and  SiO  (Goldreich  and  Kwan  1974)  are 
probably  not  important  compared  to  CO  and  IID.  Atomic  coolants  such  as 
CI,  CII,  and  01  (cf.   Penston  1970)  are  unimportant  in  dense  clouds:  due 
to  the  attenuation  of  ionizing  ultraviolet  radiation,  carbon  is  mostly 
neutral  when  n>10  cm   (Werner  1970) ,  and  atomic  carbon  and  oxygen  are 
depleted  by  chemical  reactions  in  dense  clouds  (of.  Allen  and  Robinson 
1977).  Furthermore,  the  cross  section  for  collisional  excitation  of  the 
CI  and  CII  fine -structure  levels  by  H,  is  much  lower  than  for  atomic 
hydrogen. 

The  energy  radiated  per  unit  volume  per  second  in  a  transition 
between  states  u  and  I   is 
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7 

where  nu  is  the  number  of  r.»lecules  cm   in  state  u,  E  .  is  the  energy 
difference  between  u  and  £,  A  .  is  the  Einstein  A-coefficient  (i.e. 
the  transition  probability  per  unit  time  for  spontaneous  emission) ,  and 

l-e"T*u 

v  -  *fc-  •      CA5) 

is  the  photon  escape  probability.  The  optical  depth  in  a  line  having  a 
rest  frequency  v,  and  a  thermal  Doppler  width  (Mihalas  1970) 

Av  -  (^-)  -  ,  (A6) 

is  given  by  (Penzias  1975) 

2 

e  ^Aiti     H  ,     hvn   ' 

t„„  ■  ^—      1  -  exp(-  SK)       .  (A7) 

u       Sirg^Av       L  UJ 

In  these  equations,  A  is  the  molecule's  atomic  weight,  e    and  g.   are, 
respectively,  the  statistical  weights  of  levels  u  and  I,  N.=n„R/2  is 
the  column  density  of  molecules  in  state  I  through  a  mean  path  length 

R/2,  R  being  the  cloud  radius,  c=3.0*10      cm  s       is  the  speed  of  light, 

-27 
and  h=6. 626*10        erg  s  is  Planck's  constant.     In  Eq.   (A7) ,  the  term  in 

brackets  is  the  correction  for  stimulated  emission. 

Under  steady-state  conditions,  the  relative  population  of  levels  is 

given  by 


\      CEu/gp  exp(-Eu,/kT) 


where 


1+<W<W 


Cu    ■  » V  g  if      g  <  1 


(A8) 


-  "<VV>  if      W  >  1  ^ 


65 


is  the  collisional  de-excitation  rate,  o^  being  the  collisional 
de-excitation  cross  section,  and  B  is  the  molecule's  rotational  constant. 
The  cross  sections  are  averaged  over  Maxwellian  velocity  distributions 
to  obtain  rate  constants  <av>  for  rotational  excitation.  For  simplicity, 
we  assume  that  o^  =  c£u  (Ha),  and  furthermore,  that  the  rotational  levels 
are  excited  (and  de-excited)  by  collisions  with  H-,,  the  most  abundant 
molecular  species  in  dense  clouds.  Thus,  the  term  in  brackets  in  the 
denominator  of  Eq.  (A8).  is  a  measure  of  the  deviation  from  thermal 
equilibrium:  at  higher  gas  densities,  collisional  de-excitation  dominates 
spontaneous  emission  so  that  the  levels  become  thermalized,  and  Eq.  (A8) 
reduces  to  the  Boltzmann  distribution.  When  level  populations  become 
thermalized,  the  cooling  rate  which  is  proportional  to  n  at  lower 
densities,  goes  like  n  [see  Eq.  (A4)]  .  The  cooling  efficiency  is  thus 
reduced  at  high  densities,  and  the  cooling  is  said  to  be  'coHisionally 

quenched . ' 

The  rotational  levels  of  a  diatomic  molecule  are  indexed  by  the 
rotational  angular  momentum  quantum  number  J.  The  level  J  has  energy 

Ej  =  hBJ(J+l)   ,    .  t«« 

and  the  level  degeneracy  is 

gj=2J+l       .  ^ 

The  energy  separating  levels  J  and  J-l  (corresponding  to  an  electric 
dipole  transition  governed  by  the  selection  rule  AJ=*1)  is 

EJ"EJ-1  =  ZhBJ   "  (A12) 

For  thermalized  levels,  the  population  of  the  Jth  level  is,  from  Eq.  (A8) , 
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n,  =  n0C2J+l)exp[-J(J+l)h?/kT]   ,  (A13) 

where  n  is  the  population  of  the  ground  state.  The  total  number 
density  of  a  particular  molecular  species  is 

"t  =  I    nj 


-  n  I     (2J+l)exp[-J(J+l)hE/kT]   .  (A14) 

0  J=0 

For  kT»h8   ,  as  is  the  case  for  all  molecules  considered  here, 
L'l  (2J+l)exp[-J(J+l)hB/kT]dJ 

1   °  Jo 

Using  this  result  to  eliminate  nQ  in  favor  of  n  in  Eq.  (A13) ,  we  have 

nj  »  n  H  (2J+1)  exp[-J(J+l)hS'kT]  .  (A16) 

The  cooling  rate  per  unit  volume  for  electric  dipole  transitions  is, 
after  Eq.  (A4) , 

A  "  jj  nJEJ,J-lAJ,J-l6J  •  CM7) 

with  nJ(  E     -j  and  2j  defined  by  Eqs.    (A16) ,   (A12) ,  and  (A5) ,  respec- 
tively, and 

512tt  BuJ  (A18) 

J,J  L       3hc  (2J+1) 

where  ;i  is  the  molecule's  electric  dipole  moment. 

For  lines  which  are  optically  thick,  cooling  occurs  from  the  surface 
of  a  cloud  at  a  rate  per  unit  volume  given  by 


6? 


3II.B  dv 


where  n./n  is  given  by  Eq.  (A16)  and  the  photon  spectral  energy  distri- 
bution for  i»l  is  given  by  the  Planck  function, 

Mi -^  j*n    •  '        (A20) 

For  thermal  line  broadening,  dv  is  given  by  Eq.  (A6) . 

Molecular  Hydrogen  • 

The  H,  molecule  is  a  homonuclear  diatomic  molecule  and  therefore 
has  no  permanent  electric  dipole  moment.  For  low  temperatures  (T<100  K) , 
hydrogen  molecules  are  predominantly  in  the  ground  rotational  state,  so 
that  cooling  occurs  mainly  via  the  electric  quadrapole  J=2-*0  transition 
for  para-H?.  At  such  low  temperatures,  para-ortho  collisional  conversion 
is  very  slow.  Because  of  the  long  radiative  lifetimes  of  the  excited 

levels,  the  level  populations  for  gas  densities  n>100  cm  ,  are  given 

-14 
by  the  Boltzmann  equation.  For  the  J=2-K)  transition,  E20  =  7.07x10   erg 

so  that 

n7  -  5n  e"512/T    .  (A21) 

2    o     .  . 

The  cooling  rate  by  H,,  assuming  n  =11^  =n  is  then 

L 

AH2  =  n2E20A20S2 

-  1.04xl0"23nS2e"512/T  erg  s"1  cm"3  ,  (A22) 

where  we  have  set  A  Q  =  2.947xl0"n  s"1  (Thaddeus  1972).  The  photon 
escape  probability  is  related  to  tlie  optical  depth,  as  defined  by  Eq.  (A7) , 
through  Eq.  (A5) .  For  H,, 
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t2(H2)  =  3.48*10"V*  (m/me)3  n3  ,  (A23) 

and  we  have  used  p=imJi=3m/4TrR  to  eliminate  R  in  favor  of  n  and  m,  the 
cloud  mass.  Surface  cooling  follows  from  Eq.  (A19) : 

A^  =  1.77xlO"19T%  n3  (m/m@)"3  e"512/T(l-B2)  erg  s"1  cm"3.   (A24) 

Because  of  the  low  temperatures  characterizing  dense  clouds,  KL  cooling 
is  not  expected  to  be  very  efficient. 

Hydrogen  Deuteride 

The  importance  of  HD  as  a  molecular  coolant  in  dense  clouds  was 

first  suggested  by  Dalgarno  and  Wright  (1972) .  Due  to  vibronic  inter- 

-22 
action,  HD,  unlike  H?,  has  a  permanent  dipole  moment  p=5. 85x10    esu  cm. 

For  T<100  K  only  the  first  two  excited  rotational  levels  need  be  considered. 

The  HD  cooling  rate  per  unit  volume  is  then 

V  =  jx  njEj,j-iV-A  •  CA25) 

Dalgarno  and  Wright  give  E,,  =  3.54*10"14  erg,  E.fl  =  1.77x10"   erg, 
A?,  =  2.4M0"7  s"1,  and  A  „  =  2.54x10"  s  .  The  levels  are  thermalized 
at  relatively  low  densities,  and  n  =iw,  for  Ts65  K,  as  can  be  seen  from 
Eq.  (A15) .  The  HD  cooling  rate  is  then 

A^  =  x(HD)n(1.35xl0"21B1e"129/T+1.28xl0~19B2e":585/T)  erg  s"1  cm"3  , 

(A26) 

where  xCHD)EnHn/nH  is  the  fractional  abundance  of  HD  relative  to  H.,. 

2  -5 

A  reasonable  H:D  ratio  is  20,000:1.  Accordingly,  we  take  x(HD)  =  5x10  . 

Observations  of  HD:H  by  Spitzer  et  al.    (1973),  and  DCN:HCH  by  Wilson 

et  al.    (1973)  suggest  that  x(HD)  may  be  two  orders  of  magnitude  larger. 
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However,  we  follow  the  ideas  of  Watson  (1973)  and  assume  that  these 
ratios  reflect  chemical  fractionation  rather  than  true  isotopic 
abundances.  The  photon  escape  probabilities  are  determined  by  the 
optical  depths  which  follow  from  Eq.  (A7) : 

t,(HD)  -  1.39x(HD)T  T(in/mJ3  n* 
JL  0 

x2(HD)  =  2.7Sx(HD)T"5"(m/ine)TnT  e"129/T  .  (A27) 

The  surface  cooling  rate  for  x(HD)  =  5-10"  follows  from  Eq.  (A19) ,  and 

is  given  by 

_  j_ 
J^P  =  3.46xlO"22A"5"  (m/me)  e"129/T[a-B1)+79e"256/T(l-B2)]  erg  s"1  cm"3. 

(A28) 

Carbon  Monoxide 

The  most  abundant  heavy  molecule  in  dense  interstellar  clouds  is  CO. 
Because  only  5.5  %  separates  the  first  excited  rotational  level  from  the 
ground  state,  CO  is  a  potentially  efficient  coolant  in  the  relatively 
cool  (T=10  K)  environment  of  molecular  clouds.  Indeed,  Glassgold  and 
Langer  (1973b)  have  shown,  neglecting  optical  depth  effects,  that  the 
low  temperatures  [T=10  K;  see  Zuckerman  and  Palmer  (1974)  for  references] 
typically  observed  in  dense  molecular  clouds  can  be  maintained  by  CO  cooling 
alone. 

The  CO  cooling  rate  follows  from  Eq.  (A17) : 
J. 


AC0  =  r  ' 


jj  nJ  EJ,J-1  AJ,J-16J   ' 


where  the  summation  is  discontinued  at  the  first  level  P=J_)  whore  the 
collisional  de-excitation  rate  is  less  than  the  spontaneous  transition 
rate.  The  cross  section  for  rotational  excitation  of  CO  by  H,  is 


IF     ? 

a-10        cm  (Green  and  Thaddeus  1976).  Because  of  its  low  dipole 

-19 

moment  (u=l.  12x10   esu  cm)  and  high  abundance,  the  CO  molecule 

thennalizes  at  low  densities.  Using  Eqs.  (A16) ,  CA12) ,  and  (A18) , 
Eq.  (A29)  becomes 

5    J 

ACQ  =  7.42x10"76x(CCO  ^r-    lm  Bj  J5  exp[-J(J+l)hB/kT]   erg  s"1  cm"3  . 

(A30) 

12-5 
Ke  adopt  a  fractional  abundance  x(  CO)  ■  3*10   for  the  main  isotopic 

it  n 

species,  and  x(    CO)  =  3.4x10"     for  the  less  abundant  species;  the  adopted 

12 

isotopic  ratio  is  the  terrestrial  value  89:1.  For   CO,  8=57,700  MHz, 

and  for   CO,  B=55,100  MHz.  The  optical  depth  in  a  line  arising  from  a 
transition  where  J+J-l  is,  from  Eq.  (A7)  , 

-J.     _L  A 
t,(C0)  ■  7.57xl0"6x(CO)  JT  2(m/m0)V  exp[-J(J+l)h8/kT]  • 

•[l-exp(-J(J+l)hB/kT)]   .  (A31) 

Surface  cooling  in  tlie  optically  thick  lines  occurs  at  a  rate 

JU  «  5.62xl0"82   g5"T    f  J4(2JU)(1-6T)  ' 
T*(m/mJ'  J'X 

& 


exp[-2J(J+l)hB/kT]" 
l-expl-J(J+l)hB/kT] 


erg  s"1  cm"3   .        (A32) 


Grain  Cooling 
■Cooling  by  inelastic  collisions  of  grains  with  the  gas  particles 
(mostly  molecular  hydrogen)  at  temperature  T  occurs  at  a  rate  (Dalgarno 
and  McCray  1972) 

A  =  2xlO"33n2THT-T  )0  erg  s"1  cm"3  ,  (A33) 
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where  T  is  the  temperature  of  the  grains,  and  0  is  the  energy  accom- 
modation coefficient,  and  is  a  measure  of  the  elasticity  of  the  collision; 
6  is  unity  for  a  completely  inelastic  collision,  and  becomes  zero  for 
elastic  collisions.  We  take  0=1.  Because  of  the  strong  density  dependence, 
grain  cooling  is  expected  to  dominate  molecular  cooling  at  high  densities. 
The  grain  temperatures  can  be  determined  from  the  energy  balance  of  a 
grain  (Low  and  Lynden-Bell  1976) : 

aTjQp(Tg)=nk(T-Tg)(^)0  ♦  oT^)   ,  (A34) 

where  a=S.67><10"5  erg  an   deg"  s"  is  the  Stefan-Boltzmann  constant, 
0  (T)  is  the  Planck  mean  absorption  efficiency,  n^  is  the  mass  of  the 
hydrogen  molecule,  and  T.»2.7  K  is  cosmic  blackbody  radiation  temperature. 
The  first  term  represents  the  heat  loss  due  to  radiation  of  a  dust  grain 
at  temperature  T  .  The  second  term  is  the  collisional  energy  gain  from 
the  gas,  and  the  last  term  represents  heating  from  an  isotropic  blackbody 
radiation  of  temperature  T.  .  Eq.  (A34)  assumes  that  the  cloud  is 
shielded  from  the  external  stellar  radiation  field,  and  that  there  are 
no  embedded  stars  within  the  cloud  which  might  heat  the  grains  to  a 

temperature  T  >T,  in  which  case  the  grains  heat  the  gas  (Leung  1976) . 

7  2  67 
Kellman  and  Gaustad  (1969)  give  Q  (T)  =  4.1*10"  T  '   for  0.2  micron 

ice  grains. 

The  optica]  depth  through  a  cloud  of  radius  R  for  the  thermal 

radiation  from  the  grains  is 

T  (T  )  =  0  (T  )o  n  R   .  (A35) 

g^  g     V  r  ii 

By  eliminating  R  in  favor  of  n  and  m,  the  mass  of  a  uniform,  spherical 
cloud,  and  by  taking  a  n    =  6*10"22n,  and  Q  (T  =10)=10"4,  the  optical  depth 

*»    4>  it       h 
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T  (T  =10)  =  10"7(m/rae)Tn3  .  (A36) 

For  stellar-mass  clouds,  » 361  for  densities  nSlO   cm  ,  which  is  in 
good  agreement  with  the  results  of  Hattori  et  al.    (1969).  Since  densities 
encountered  during  the  initial  collapse  stages  are  much  less  than  this, 
the  thermal  radiation  from  grains  is  assumed  to  pass  freely  out  a 
contracting  cloud. 

Cloud  Temperature 

Figure  Al  illustrates  the  temperature  and  density  dependence  of  the 
various  heating  and  cooling  rates  described  in  this  Appendix.  The  total 
molecular  cooling  rate  7^=^  +AHD+AC0  •  As  expected  grain  cooling  com- 
pletely dominates  for  n>106  cm"3.  Even  if  gas-grain  collisions  are  only 
weakly  inelastic  (0=0.01  say),  grain  cooling  will  still  dominate  at  high 
densities.  Grain  temperatures  ranged  from  T  =6  K  for  WO  K  up  to 
T  =35  K  for  T=100  K.  Gas  temperatures  much  in  excess  of  20  K  are  rarely 
encountered  in  dense  clouds;  temperatures  up  to  100  K  are  considered 
merely  to  illustrate  the  rapid  rise  in  the  molecular  cooling  rates  at 
these  relatively  high  temperatures. 

Surface  cooling  for  all  molecules  never  amounts  to  more  than  10%  of 
the  total  molecular  cooling  rate,  and  cooling  by  molecular  hydrogen  becomes 
noticeable  (-101  of  the  total  molecular  cooling)  only  for  the  highest 
temperatures.  Carbon  monoxide  dominates  the  cooling  at  T=10  K  for  all 
densities.  However,  except  for  n=10  cm   where  KwT^ro'  h'Jdr°9en  deuteride 
is  by  fax  the  dominant  molecular  coolant  for  T  2  20  K,    its  cooling  rate 
becoming  three  orders  of  magnitude  greater  than  CO  cooling  for  the 
highest  temperature. 

.  .   It  appears  that  cloud  temperatures  equivalent  to  those  typically 
observed  in  dark  clouds  (T=10  K;  see  Zuckerman  and  Palmer  (1974)  for 
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references)  are  certainly  attainable  at  low  gas  densities,  and  should 

prevail  for  higher  densities  as  well,  provided  grain  cooling  operates 

with  at  least  a  one-percent  efficiency.  The  rapid  rise  of  the  molecular 

cooling  rate  at  higher  temperatures  indicates  that  a  cloud  will  collapse 

approximately  isothermally  near  some  equilibrium  (r=A)  temperature 

T  =10-20  K  at  least  for  densities  n<10  cm-  .  Molecular  cooling  becomes 
e 

even  more  efficient  if  large  velocity  gradients  (e.g.  due  to  cloud  collapse) 
develop  in  a  contracting  cloud  (Goldreich  and  Kwan  1974;  de  Jong  et  al. 
1975).  Furthermore,  as  is  shown  in  Appendix  B,  the  cosmic-ray  ionization 
rate  (and  therefore,  the  heating  rate)  adopted  here  is  probably  greatly 
overestimated  for  densities  n>10  cm   .     --  — , 


APPENDIX  B 
FRACTIONAL  IONIZATION  IN  DENSE  MAGNETIC  CLOUDS 

The  degree  of  ionization  in  a  dense  interstellar  cloud  controls  the 
coupling  of  the  cloud  to  the  galactic  magnetic  field.  When  the  fractional 
ionization,  n  /n,  is  high,  the  time  scale  for  ambipolar  diffusion  is 
long. 

For  an  ionizing  flux  C ,  the  electron  density  ng  in  a  cloud  having 

a  gas  density  n  varies  with  time  according  to  (Oppenheimer  and  Dalgarno 

1974) 

dn  +        +        + 

-J^  =  £n-a  n(m  )n  -arn(M  )ng-a  n(M  )n   ,  (Bl) 

where  n(m+)  and  n(M+)  denote  the  density  of  molecular  and  heavy  metal 


ions,  respectively,  and  a,  «  ,  and  a     are,  respectively,  the  rate 

+ 
coefficients  for  dissociative  recombination  of  molecular  ions  m  , 

radiative  recombination  of  heavy  metal  ions  M  ,  and  recombination  on 
grains.  For  simplicity,  it  is  assumed  that  all  rate  coefficients  are 
independent  of  ion  type,  It  is  easily  shown  that  for  densities 
n>10  cm"  ,  recombination  on  grain  surfaces  dominates  radiative  recom- 
bination so  that  the  third  term  on  the  right -hand -side  of  Eq.  (Bl)  can  be 
neglected.  In  steady-state  equilibrium,  Eq.  (Bl)  then  gives 

&!-an(m+)n 
n(M  )  = .  (B2) 


The  electron  density  is 

ne  =  n(m+)  +  n(M+)   .  (B3) 
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The  equilibrium  rate  equation  for  the  molecular  ion  density  gives 

(Oppenheimer  and  Dalgamo  1974) 


n(m  )  ■ 


Sn 
un  +gn(M) 


(B4) 


where  n(M)  is  the  total  density  of  heavy  neutral  atoms  (e.g.  Mg,  Ca,  Na, 
and  Fe)  that  undergo  charge  transfer  with  molecular  ions,  and  8  is  the 
rate  coefficient  for  the  charge-transfer  process.  We  neglect  associative 
ionization  (Oppenheimer  and  Dalgarno  1977)  which  is  probably  an  insig- 
nificant source  of  electrons  in  cold  interstellar  clouds. 

Putting  Eqs.  (B2)  and  (B4)  into  Eq.  (B3)  gives  a  quadratic  equation 


for  n  which  has  the  solution 
e 


_  Bn(M)  +  1 


Bn(M) 
a 


4£ 
a 


Bn(M) 


-  "g 


1/2 


(B5) 


2a     2 

Following  Oppenheimer  and  Dalgamo  (1974) ,  we  define  a  depletion  factor 
6  by  the  expression 

(B6) 


n(M) 


4*lCf56n 


Observations  with  the  Copernicus  satellite  (Morton  et  al .   1973;  Morton 

1974,  1975)  suggest  an  average  depletion  factor  for  heavy  metals,  6=0.1. 

Oppenheimer  and  Dalgamo  give  6=10   cm  s  ,  a=10   cm  s  ,  and 

a  =10-17  cm3  s"1. 
r 


Eq.   (B6)  then  becomes 


ne=2xl0-9n     ^J^l/2.! 


so  that  for  densities  n»10    5, 


V«ft 


(B7) 


(B8) 
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Ions  in  dense  molecular  clouds  are  supplied  primarily  by  the 

40 
ionization  of  molecular  hydrogen  by  cosmic  rays  and   K  radioactivity; 

ultraviolet  radiation  (Werner  1970)  and  X-rays  (Nakano  and  Tademaru 

3   -3 

1972)  are  screened  in  the  peripheral  regions  of  dense  (n>10  cm  ) 

clouds.  Cameron  (1962)  has  estimated  the  ionization  rate  by  the  6-decay 

40 
of   K  radioactive  nuclei  to  be  . 

5k  '   1.4*10"21  s"1  .  <B9) 

Nakano  and  Tademaru  (1972)  have  calculated  the  ionization  rate  of  atomic 
hydrogen  by  cosmic  rays.  Adjusting  their  result  for  ionization  of 
molecular  hydrogen  by  taking  into  account  the  difference  in  ionization 

cross  sections,  a„  =1.65a„  (Bates  and  Griffing  1953),  the  cosmic -ray 

2  * 

ionization  rate  in  a  cloud  having  a  mass  m  is  given  by 

KCR  =  10"17  exp[-1.54x]0"7(m/ms)~nT]  .  (BIO) 

The  exponential  term  reflects  the  attenuation  of  cosmic  rays  due  to 
their  interaction  with  matter  as  they  propagate  through  a  cloud.  The 
total  ionization  rate 


I  -  5CR  ♦  5k  •  '  (Bll) 

By  comparing  Eqs.  (B9)  and  (BIO),  one  can  see  that  for  stellar-mass  clouds 

10   12   -3 
having  densities  njlO  -10   cm  ,  cosmic  rays  are  effectively  screened  and 

40 
ions  are  produced  primarily  by  in  situ      K  nuclei. 


Brown  and  Marcher  (1977)  have  shown  that  ionization  of  H  and  R,  in  ' 
dense  clouds  may  be  enhanced  by  energetic  secondary  electrons  produced  by 
knock-on  collisions,  neutron-decay  reactions,  and  pion-decay  reactions, 
l**^****,  following  the  interaction  of  fast  (primary)  cosmic  rays  with  the 
material  in  a  cloud.  This  effect,  which  may  be  significant  only  when  low- 
energy   cosmic  rays  are  excluded  from  dense  clouds,  is  difficult  to  estimate 
quantitatively,  and  is  therefore  neglected  here. 
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Eq.  (BIO)  may  greatly  overestimate  the  cosmic-ray  ionization  rate 
ta  dense  magnetic  clouds.  If  the  magnetic  field  within  a  contracting 
cloud  becomes  tangled  (e.g.  by  turbulence),  cosmic  rays,  constrained  to 
move  along  the  magnetic  lines  of  force,  random  walk  through  the  cloud 
and  must  traverse  more  matter  to  reach  the  central  regions  of  the  cloud. 

Nakano  and  Tademaru  (1972)  have  estimated  the  importance  of  this  effect, 

S   -3 

and  have  concluded  that  for  densities  n>10'  cm  ,  cosmic  rays  are 

i 

effectively  screened,  i.e.,  t_<t  so  that  S^Sj,  . 

When  cosmic  rays  stream  along  magnetic  fields  in  a  collisionless 
plasma  faster  than  the  Alfven  velocity,  they  generate  hydromagnetic 'waves 
which  in  turn  scatter  the  cosmic  rays  (Wentzel  1974,  and  references  cited 
therein) .  Indeed,  Skilling  and  Strong  (197C)  have  shorn  that  the  incoming 
cosmic-ray  flux  in  a  dense  cloud  may  be  substantially  reduced  by  this 
mechanism.  However,  the  damping  effect  of  ion-neutral  collisions  may 
inhibit  the  generation  of  such  waves  by  the  cosmic  rays  themselves 
(Kulsrud  and  Pearce  1969) . 

As  pointed  out  by  Nakano  and  Tademaru  (1972),  cosmic  rays,  although 
not  unstable  to  the  generation  of  hydromagnetic  waves,  are  strongly 
influenced  by  the  presence  of  such  waves  generated  by  other  mechanisms 
(e.g.  by  magnetic  braking  during  cloud  collapse) .  Cosmic  rays  are 
scattered  by  these  waves  and  can  not  freely  stream  along  the  open  magnetic 
field  lines.  Nakano  and  Tademaru  have  shown  that  effective  screening  of 
the  cosmic  rays  occurs  if 

<6B2>  >  iO   ,  (B12) 


where  <6B  >/8ir  is  the  energy  density  of  the  hydromagnetic  waves-  having 
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amplitude  SB,  B  is  the  strength  of  the  cloud's  large-scale  uniform  field, 

and  p=3m/4nR  is  the  mass  density  of  the  cloud  material.  By  eliminating 

-2 
p  in  favor  of  m  and  R,  and  by  defining  n=R/R  with  B=B  n   (see  Section  III 

and  IV  of  the  text  for  a  discussion  of  the  rate  of  growth  of  the  magnetic 

field  in  a  contracting  protostar) ,  with  subscripts  denoting  initial  values, 

this  expression  becomes 

,   8*10~24R  B 
<6BS> °-°  .  (B13) 

(m/m  )n 
© 

12 

From  Table  1  in  the  main  text,  we  see  that  R  B  *  constant  =  9x10      cm  gauss 

for  all  cloud  masses.     Thus,  if 

<6B2>>^i5 ,  (B14) 

(m/m0)n 

cosmic  rays  are  effectively  screened. 

The  energy  lost  from  a  magnetically  braked,  collapsing  cloud  is  just 
the  difference  between  the  cloud's  rotational  kinetic  energy  given  by 
angular  momentum  conservation,  and  its  magnetically-braked  rotational 
energy.  For  a  cloud  which  is  constrained  to  co-rotate  with  its  sur- 
roundings, this  energy  is  given  by 

E,  _  o  O.anAV2  •  (El 5) 

lost       0  0 

Assuming  that  this  energy  is  carried  away  by  the  hydromagnetic  waves  which 
are  generated  by  the  braking  process,  we  can  write 

<6B2>=2xlO"11  n"5  ,  (B16) 

2  3 
where  we  have  taken  initial  values  from  Table  1,  noting  that  to R  /(m/m  ) 

25   3-2 

=  constant = 10   cm  s  .  Actually,  depending  on  the  mass  of  the  cloud, 

7 
<6B  >  may  be  an  order  of  magnitude  smaller  or  larger,  so  that  Eq.  (B16) 

represents  an  approximate  mean  value.  A  comparison  of  Eqs.  (B14)  and 
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(B16)  shows  that  cosmic  rays  are  effectively  excluded  from  a  contracting 
cloud  very  soon  after  the  collapse  begins. 

It  is  for  these  reasons  that  we  believe  the  ionization  in  dense 
oagnotic  clouds  is  determined  by   the       K  ionization  rate,    I.e.,  C*"?^  • 
This  being  the  case,  the  magnetic  flux  linking  a  contracting  cloud  to 
its  surroundings,  uncouples  at  relatively  low  gas  densities.  This  may 
explain  the  absence  of  large  magnetic  fields  in  some  dust  clouds 
(Crutcher  at  al .   1975) . 

Even  if  the  cosmic  rays  are  not  magnetically  scattered,  there  will 
be  a  reduction  in  the  flux  of  cosmic  rays  in  a  magnetic  cloud.  The 
magnetic  field  lines  in  the  neighborhood  of  a  contracting  cloud  diverge 
outward  from  the  cloud  so. that  charged  particles  streaming  along  the  field 
lines  into  the  cloud  will  be  reflected  by  the  'magnetic  mirror  effect.'  . 
Fermi  (1949,  1954)  proposed  that  such  a  magnetic  reflection  mechanism 
might  explain  the  origin  of  the  galactic  cosmic  rays.  For  slow  variations 
of  the  magnetic  field  in  time  and  space,  the  diamagnetic  moment  of  a 
charged  particle  is  an  adiabatic  invariant.  Let  8  be  the  angle  between 
the  direction  of  the  line  of  force  and  the  direction  of  motion  of  the 
spiraling  particle,  viz.,   the  pitch  angle.  Assuming  an  isotropic  distri- 
bution of  particle  velocities  in  a  region  where  the  field  strength  is  B, 
one  can  easily  show  (cf.   Spitzer  1962)  that  the  velocities  must  fall 
within  a  solid  angle  defined  by  the  pitch  angle  such  that 

8  =  sin^Wn/2  CB17) 

when  the  field  strength  increases  to  Bg.  In  our  case,  B  is  taken  to  be 
the  (uniform)  galactic  field  far  from  a  contracting  cloud,  and  Bs  is  the 
field  strength  at  the  cloud  surface.  Assuming  that  the  particle  density 
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in  a  given  region  is  proportional  to  the  size  of  the  solid  angle  given 
by  Eq.  (B17),  it  follows  (cf .   Kaplan  and  Pikelner  1970)  that  the 
fractional  decrease  in  cosmic-ray  flux  is 

f  =  1-cosB   .  (B18) 

The  cosmicray  ionization  rate  is  then  given  by 

t 

where  ?rR  is  given  by  Eq.  (BIO).  Eq.  (B19)  provides  a  workable  upper- 
limit  to  the  ionization  rate  in  dense  magnetic  clouds,  the  lower-limit 
determined  by  C,  being  the  most  likely  for  reasons  already  discussed. 
Eq.  (B17)  is  valid  only  if  the  Larmor  radius 

*.*3*S£  [C*i*-l]\-  CB20) 

of  a  cosmic  ray  having  a  kinetic  energy  E,  (MeV) ,  and  moving  in  a  magnetic 
field  B  (microgauss) ,  is  much  less  than  the  radius  of  the  cloud.  One 
finds  rL(2MeV>3.4xi010  cm  and  r,(10  GeV)=6.1xl012  cm,  so  that  in  fact, 
r,«R.  Also,  the  quantity  (g-  |jp|)   must  be  significantly  greater  than 
the  Larmor  period.  It  is  easily  demonstrated  that  this  condition  is 
equivalent  to  the  dynamic  time  scale  being  much  larger  than  the  Larmor 

period,  or 

1.4(E,+938) 

tf  »  *g= s   .  (B21) 

V 

Since  collapse  times  are  typically  10  -10   s,  this  condition  is  easily 

satisfied.  Finally,  Eq.  (B17)  neglects  collisions  among  the  cosmic  rays 

which  have  the  effect  of  randomizing  the  pitch  angle  6.  We  require 

that  the  time  between  collisions  be  much  greater  than  the  Larmor  period, 
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or  equivalently,  that  the  density  of  material  (primarily  molecular 
hydrogen)  in  a  cloud  satisfies: 


14B 
U 


5.1x10  "*.      fL 

[(1+  gsfr-l]2  cm-5   ,  (B22) 


where  a  is  the  collision  cross  section.  Low-energy  cosmic  rays  (E,=2) 
interact  with  the  cloud  material  primarily  by  ionizing  molecular 

hydrogen.  Glassgold  and  Langer  (1973a)  give  o(2MeV)=1.89xlO~17  cm2. 

3 
High-energy  cosmic  rays  (E,>10  )  interact  mainly  by  p-p  scattering  and 

pion-production  reactions  with  a  cross  section  o(10GeV)=2xlCf   cm  (see 

Nakano  and  Tademaru  1972  for  references).  Thus,  from  Eq.  (B22) i 

n(2MeV)«4x!0  B  and  n(10GeV)«2xlO4  B  .  Since  initial  values  for  B 

U  u  U 

3  2/3 

range  from  3  up  to  10  (see  Table  1)  and  B  increases  as  n  '  during 

gravitational  collapse  (see  Section  IV  of  the  text) ,  these  conditions  are 

easily  satisfied. 


APPENDIX  C 

MOMENT  OF  INERTIA  FOR  DIFFERENTIALLY 
ROTATING  MAIN-SEQUENCE  STARS 

The  angular  momentum  of  a  rotating  spherical  body  having  a  radius  R, 

a  radial  density  distribution  p(r),  and  an  angular  velocity  field  w(r,6) 

is 

rR  p/2     r    %_4_ 

'0  'o 


j  =  *L  u(i,e)p(T)r  sin0d8dr     .  (CI) 

3  Jn  Jn 


Assuming  u(r,e)=w(r) ,  this  becomes 

f 


Stt 
3 


rR  a 

u>(r)P(r)rch-     .  (C2) 

Jn 


The  actual  distribution  of  mass  throughout  chemically  homogeneous  stars 
which  are  not  completely  convective  often  approximates  that  in  the 
'standard  model1  of  Eddington  (1926),  which  is  just  a  polytrope  of  n-3. 
For  a  rigidly-rotating  star,  tu(r)=constant=u)R,  the  angular  velocity  at 
the  surface  (r=R) .  The  integral  in  Eq.  (C2)  is  then  most  easily- 
evaluated  with  the  aid  of  the  Emden  solutions  for  an  n=3  polytrope 

(see,  for  example,  Chapter  23  of  Cox  and  Giuli,  1968).  Writing  the 

2 
angular  momentum  in  terms  of  the  moment  of  inertia  I=<MR  ,  where  k   is  the 

gyration  constant,  we  have 

J=I<"R 

=  tcMR2^   ,  (C3) 


85 


so  that  for  the  case  of  rigid  rotation,  we  find 

K  =  0.08  .  (C4) 

The  pre-main-sequence  evolutionary  models  of  Bodenheimer  and 
Ostriker  (1970)  for  rapidly  rotating  massive  stars  predict  a  marked 
differential  rotation,  with  the  central  angular  velocity  uc  being  a 
factor  of  ten  greater  than  uR.  Their  differentially  rotating  confi- 
gurations are  stable,  according  to  the  criterion  developed  by  Goldreich 
and  Schubert  (1967) .  From  Figure  5  of  Bodenheimer  and  Ostriker  (1970) 
we  approximate  the  angular  velocity  as 

u>(x)  =  10  wRe  '     , 

where  xSr/R,  and 

a  =  1.4  for  0  *  x  -  0.3 
-  1.1  for  0.3  -  x  *  1.0  . 

The  mass  distribution  for  a  polytrope  of  n=3  can  be  approximated  as 

,  ,    -  -bx 

PW  ■  P0P  e 


(CS) 


where 

b  =  20  and  c  ■  2    for  0  -  x  -  0.3 
b  =  11  and  c  =  1.5   for  0.3  -  x  -  1.0 


(Co) 


Here,  p  =  3m/4irR3  and  p  (n=3)  -  54.18  g  an"   .     Substituting  the  above 
expressions  for  u(x)  and  p(x)   into  Eq.    (C2)  ,  and  evaluating  the  integral 
using  Simpson's  Rule,  we  find 

k  =  0.28     .  (C7) 
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Die  gyration  constant  (and  therefore,  the  angular  momentum)  of  the 
differentially  rotating  configuration  is  thus  three  times  that  of  a 
rigidly-rotating  body. 
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